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INTRODUCTION 


1.  SOLID  STATE  DEVICE  RESEARCH 

Photovoltaic  detectors  with  long-wavelength  cutoffs  up  to  20  g  at  77°K  and  up  to  30  g  at  I2°K  have  been 

produced  from  annealed  Bridgman-grown  and  vapor-grown  Pb,  ,Sn  Te  crystals  in  the  0  <  x  <  0.  25com- 

1  x  x  *  |Q  | /9 

position  range.  Responsivities  as  high  as  190  V/W,  detectivities  ot  nearly  10  cm/W  S'*c  ,  and  ex¬ 
ternal  quantum  efficiencies  up  to  37  percent  have  been  obtained  at  wavelengths  up  to  12  g  in  diodes  op¬ 
erated  at  77°K.  Detectivities  at  12  K  were  generally  about  an  order  of  magnitude  higher  than  at  77°K. 

Photoconductivity  at  wavelengths  up  to  15  g  has  been  observed  at  77°K  in  Bridgman-grown  and  subse¬ 
quently  annealed  crystals  of  Pb,  Sn  Te  which  had  carrier  concentrations  as  low  as  2  10^  cm  ^ 

x  9  89  1/9 

and  77°K  mobilities  of  about  3  ■  104  cm“/V  sec.  Detectivity  values  between  10  and  10  cm/W  sec 

at  77°K  and  10^  cm/W  sec  at  4.2°K  were  measured.  Photoconduct ive  lifetimes  of  about  10  ^  and 

10  ^  sec  were  obtained  at  the  two  respective  temperatures. 

Diode  lasers  with  low  thresholds  and  emission  wavelengths  up  to  28.  1  g  have  been  fabricated  from 
vapor-grown  single  crystals  of  Pbj^Sn^Te  containing  Bi  as  the  dominant  donor  impurity.  For  x  ^.0.2, 
77°K  threshold  current  densities  were  less  than  1500  A  cm  With  pulse  bias,  a  peak  power  output 
of  0.1  W  and  an  external  quantum  efficiency  rj  *  of  0.08  were  observed  in  a  Pb^  ^SnQ  ^ ^ e  diode  at 

I2°K,  and  a  peak  power  of  0.03  W  with  77  0.01  was  measured  at  77°K.  With  CW  operation  at  12°K, 

-3  X 

5  x  10  W  with  77  0.03  was  obtained. 

'  x 

Extrinsic  far-infrared  photoconductivity  has  been  observed  at  4.2°K  in  high-purity  n-tvpe  epitaxial 

3 

layers  of  GaAs  grown  on  Cr-doped  semi-insulating  GaAs  substrates.  Responsivities  as  high  as  5  x  10 
4 

and  4  x  i()  V/W  have  been  measured  at  195  and  337  g,  respectively.  At  902  g,  the  responsivity  dropped 
9 

to  9  x  K)“  V/W.  Measurements  of  the  noise  in  the  detection  system  at  300  Hz  in  a  1-llz  bandwidth 
yield  an  NEP  of  1.2*  l(f  ‘  1  W  at  195  n,  1.4  -  10'12  W  at  337  n,  and  6  x  l(f 11  W  at  902  n.  The  time 
constant  of  the  detector  has  been  determined  to  be  shorter  than  l  gsec  using  a  Ge  avalanche  modulator 
to  chop  the  incident  radiation.  A  time  constant  of  about  5  nsec  was  measured  using  impact  impurity 
ionization  in  the  GaAs. 

We  have  found  that  proton  radiation  damage  can  be  used  to  convert  both  p-  and  n-type  GaAs  into  high- 
resistivity  material,  and  have  used  this  technique  to  isolate  p-n  junctions  on  a  diffused  GaAs  substrate 
and  also  to  prevent  edge  breakdown  in  Au-GaAs  Schottky  barriers.  This  technique  of  creating  high- 
resistivity  material  should  be  useful  for  the  fabrication  of  structures  in  which  thin  high-resistivity 
layers  are  required. 

The  photovoltaic  response  oflnSb-MOS  detectors  to  modulated  3.9-g  radiation  from  anlnAs  diode  emit¬ 
ter  has  been  measured  at  77°K  following  irradiation  with  short  wavelength  light  in  the  energy  range  be¬ 
tween  1  and  5  eV.  Effects  are  observed  due  to  trapping  and  detrapping  of  electronic  charge  at  oxide 
states  and/or  at  oxide-InSb  interface  states;  however,  the  precise  model  for  this  trapping  has  not  yet 
been  determined. 
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II.  OPTICAL  TECHNIQUES  AND  DEVICES 

The  100-W  C02  laser  amplifier  has  been  operated  in  the  sealed-off  mode  and  yields  an  output  of  88  W 
compared  with  1 14  W  using  a  flowing  gas  mixture.  The  output  beam  pattern  is  a  close  reproduction  of 
the  input  TEMqq  mode. 

The  1.5-meter  stable  CC>2  oscillator  has  been  used  to  drive  the  above  amplifier.  An  output  power  of  10 
to  15  W  is  obtained,  and  the  frequency  stability  should  be  comparable  to  that  of  previous  smaller  ver¬ 
sions  once  the  output  window  heating  is  eliminated  by  the  use  of  lower-loss  materials. 

An  attempt  has  been  made  to  observe  thermal  blooming  of  a  10.6-g  C02  laser  beam  in  an  absorbing 
gas.  Substantial  blooming  does  occur,  but  difficulty  has  been  encountered  in  making  quantitative  meas¬ 
urements  due  to  the  asymmetry  caused  by  convection  in  our  horizontal  gas  cell  and  to  the  lack  of  a  good 
image  detector  at  the  10 -g  wavelength. 

Optical  heterodyning  at  10.  6  gm  between  a  current-tunable,  single-mode  Pb^  ggSn^  ^Te  diode  laser 
and  a  stable,  single-frequency  C02  gas  laser  has  produced  a  controllable  beat  frequency  in  a  Ge:Cu  de¬ 
tector  from  under  50  to  1300  MHz.  From  these  measurements,  we  have  determined  the  cutoff  frequency 
of  a  Ge:Cu  detector  considered  for  application  in  the  laser  radar  system,  and  have  obtained  a  direct 
indication  of  the  linewidth  of  emission  from  the  semiconductor  p-n  junction  laser. 

III.  MATERIALS  RESEARCH 

A  theoretical  analysis  of  the  forced  convection  method  of  vapor  crystal  growth  has  shown  that  the  growth 
velocity  is  proportional  to  the  first  power  of  the  available  concentration  of  condensing  vapor  and  to  the 
square  root  of  the  carrier  gas  velocity.  These  results  are  supported  by  experimental  data  reported 
previously  on  the  growth  of  iodine  and  camphor  crystals.  The  condition  required  for  stability  of  the 
growth  interface  has  been  formulated  in  terms  of  a  relationship  between  the  temperatures  at  the  inter¬ 
face,  at  the  edge  of  the  diffusion  boundary  layer  adjacent  to  the  interface,  and  at  the  solid  source. 

Optical  transmission  measurements  have  been  made  at  room  temperature  on  single  crystals  of 
Cdj_xMg^Te  alloys  with  0  <  x  ^  0.46.  The  energy  gap  evaluated  from  the  data  increases  linearly  with 
increasing  x,  from  1.50  eV  for  CdTe  to  2.23  eV  for  x  =  0.  46. 

Photoluminescence  spectra  due  to  oxygen  in  ZnTe,  Se  ,  ZnTe,  S  and  Zn,  Cd  Te  alloys  have  been 

1  X  X  1  -X  X  1  -X  X 

measured  at  4.2°K.  As  expected  for  the  isoelectronic  trap  model  of  the  oxygen  center,  the  trapping 
energy  decreases  with  increasing  ZnSe  and  ZnS  content,  but  is  essentially  independent  of  CdTe  content. 

The  pressure-composition  phase  diagram  for  the  Ba^S^RuO^  system  at  1000°C  has  been  determined 
by  x-ray  diffraction  measurements  on  powdered  samples.  Four  different  structures,  with  different 
proportions  of  cubic  and  hexagonal  close  packing,  are  observed.  The  addition  of  strontium  and  the  ap¬ 
plication  of  pressure  both  favor  structures  with  more  cubic  and  less  hexagonal  close  packing. 

A  method  has  been  developed  for  evaluating  absolute  atomic  form  factors  in  solids  by  analysis  of  inte¬ 
grated  intensity  data  obtained  in  x-ray  diffraction  measurements.  Application  of  this  method  to  data  for 
nickel  suggests  a  transfer  of  0.3  electron  from  the  outer  levels  (4s)  to  the  first  inner  levels  (3d),  but 
this  should  be  considered  as  a  qualitative  result. 
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IV.  PHYSICS  OF  SOLIDS 

The  optical  investigation  of  the  band  structure  of  nickel  has  now  been  extended  to  nickel-copper  alloys. 
Reflection  measurements,  from  0.2  to  11  eV,  in  several  alloys  seem  to  indicate  that  below  4  eV  the 
alloys  behave  according  to  a  "localized  states"  or  to  a  "minimum  polarity"  model. 

The  study  of  the  effect  of  stress  on  the  impurity  spectrum  of  bismuth-doped  silicon  lias  now  been  com¬ 
pleted.  All  measured  relative  intensities  and  transition  energies  as  a  function  of  stress  are  in  agree¬ 
ment  with  predictions  of  a  strain  Hamiltonian  derived  from  effective  mass  and  deformation  potential 
theory. 

Two  sets  of  interband  transitions,  analogous  to  those  in  HgTe,  have  now  been  observed  in  the  magneto¬ 
reflection  of  gray  tin.  One  of  the  lowest  energy  transitions  does  not  fit  into  the  Landau. level  scheme; 
it  is  suggested  that  this  line  may  arise  from  transitions  between  impurity  states  associated  with  the 
Landau  levels . 

Investigation  of  the  "extra"  transitions  in  the  interband  magnetoreflection  and  magneto-electroreflection 
of  InSb,  caused  by  the  warping  and  linear-k  splitting  of  the  valence  band,  has  been  concluded.  The 
warping  obtained  is  in  good  agreement  with  the  cyclotron  resonance  results  of  Bagguley,  etal.;  the 
linear-k  splitting  is  one-third  the  value  deducedfrom  cyclotron  resonance  by  Robinson,  but  is  in  agree¬ 
ment  with  the  theoretical  estimate  of  Kane. 

A  technique  has  been  developed  for  determining  the  Debye  temperature,  0^,  of  crystals  by  means  of 
thermal  modulation  of  x-ray  diffraction.  This  new  method,  which  has  been  tested  on  NaCl,  gives  re¬ 
sults  in  good  agreement  with  other  techniques,  and,  aside  from  holding  promise  of  surpassing  the  ac¬ 
curacy  of  existing  DC  methods,  has  the  additional  advantage  of  yielding  0^  from  measurements  at  a 
single  temperature. 

Substantial  amplification,  greater  than  25db/cm,  of  9-GHz  longitudinal  ultrasonic  waves  has  been  ob¬ 
tained  in  n-InSb  on  application  of  a  pulsed  electric  field.  Theoretical  calculations  of  this  effect  are  in 
reasonable  qualitative  agreement  with  the  experiment. 

Measurements  of  the  magnetic  properties  of  the  high-pressure  (orthorhombic)  phase  of  MnAs  have  been 
extended  to  11  kbars.  The  variation  with  temperature  of  the  magnetic  moment  has  been  studied  as  a 
function  of  magnetic  field  and  throughout  the  pressure  range. 

A  study  of  the  properties  of  the  one-dimensional,  two-sublattice  model  of  a  spiral  magnetic  structure 
has  been  completed.  For  a  range  of  parameters,  a  ferrimagnetic  conical  spiral  structure  is  the  ground 
state.  The  spin  wave  spectrum  for  this  case  has  been  investigated. 

The  well-known  spherical  model  of  a  ferromagnet  has  been  shown  to  be  thermodynamically  identical  to 
a  system  of  infinite-dimensional  spins  interacting  through  a  Heisenberg  Hamiltonian. 

It  is  known  that  for  the  Ising  model,  the  onset  of  magnetic  order  is  characterized  by  the  collapse  of  the 
eigenvalue  spectrum  of  a  certain  linear  operator.  It  has  been  demonstrated  that  the  same  mathemati¬ 
cal  phenomenon  characterizes  the  onset  of  ordering  in  a  set  of  y-dimensional  spins  interacting  via  a 
Heisenberg  Hamiltonian. 
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A  new  form  of  temperature-dependent  Hartree-Fock  theory  is  being  investigated,  based  on  a  variational 
principle  for  the  free  energy.  The  one-electron  orbitals  used  to  construct  single  determinant  wave 
functions  must  satisfy  an  equation  which  is  more  general  than  the  usual  Hartree-Fock  equations.  This 
new  formulation  eliminates  some  of  the  paradoxical  results  of  ordinary  temperature-dependent  Hartree- 
Fock  theory. 

Using  high-resolution  Brillouin  scattering  techniques,  the  velocity  and  attenuation  of  28-GHz  longitudi¬ 
nal  phonons  propagating  along  the  x-axis  ofcv-quartz  have  been  measured  in  the  temperature  range  300° 
to  600°  K.  Although  the  velocity  variation  with  temperature  agrees  with  the  earlier  ultrasonic  meas¬ 
urements,  the  hypersonic  attenuation  agrees  with  theoretical  extrapolations  of  ultrasonic  data  only  in 
order  of  magnitude. 

For  the  first  time,  light  scattering  has  been  observed  from  screened  single-particle  electron  and  also 
hole  excitations.  Scattering  from  electron  excitations  was  found  in  GaAs,  InP,  AlSb,  and  CdTe,  and 
from  hole  excitations  in  p-type  GaAs  doped  with  cadmium,  zinc,  and  magnesium. 

The  study  of  thermal  defocusing  of  a  laser  beam  has  now  been  extended  to  steady  state  (and  approach 
to  steady  state)  conditions.  Competing  conduction  and  convection  processes  have  been  investigated  in  a 
number  of  new  experiments. 

Methods  of  achieving  single-mode  operation  of  a  high-power  pulsed  ruby  laser  have  been  explored.  In 
the  best  arrangement,  a  single  transverse  mode  and  less  than  six  longitudinal  modes  are  obtained. 
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I.  SOLID  STATE  DEVICE  RESEARCH 


A.  PHOTOVOLTAIC  DETECTORS  OF  Pb^Sryre 

Photovoltaic  detectors  with  significantly  improved  sensitivity  have  been  fabricated  from 
vapor-grown  Pb^  ^Sn^Te  with  compositions  between  x  =  0.17  and  0.20,  for  which  the  cutoff  wave¬ 
length  at  77°K  falls  in  the  8-  to  14-p  atmospheric  window.  The  detectivity  has  been  improved 
by  increasing  both  the  quantum  efficiency  and  the  zero-bias  impedance  of  the  diodes.  Most 
diodes  were  produced  by  a  combined  operation  of  vapor-growth,  annealing,  and  diffusion  which 
was  performed  in  an  evacuated  ampoule  by  varying  the  temperature.  The  samples  were  then 
cleaved,  contacts  were  made  by  gold-plating,  and  leads  were  attached  with  a  Ga-In  alloy.  The 
junctions  were  5  to  40  p  deep.  In  most  cases,  the  deeper  junctions  had  higher  zero-bias  imped¬ 
ances.  This  seems  to  correlate  with  the  absence  of  metallic  inclusions  from  the  bulk  of  vapor- 
grown  crystals,  as  opposed  to  regions  near  the  surface  where  high  densities  of  inclusions  have 
been  detected  by  etching.  In  order  to  remove  damaged  layers  from  the  cleaved  faces  and  to 
reduce  the  junction  depth,  the  diodes  were  stream-etched  with  an  electrolytic  solution  consisting 
of  20 g  KOH,  45  ml  H^O,  35  ml  glycerin,  and  20  ml  of  ethanol. 


TABLE 

1-1 

PROPERTIES  OF  Pb. 

1 

Sn  Te  DETECTORS 

K  X 

Area 

(mm^) 

Temperature 

(°K) 

Cutoff 

Wavelength 

(m) 

Incremental 

Resistance 

(ohms) 

Efficiency 

Peak 

Responsivity 

(V/W) 

Peak 

Detectivity 

(cm/W  sec 

0.1 

77 

9.5 

180 

0.  14 

190 

5. 3  X  109 

0.25 

77 

11 

16 

0.17 

27 

9 

4X  10 

0.2 

77 

10 

10 

0.33 

26 

3. 5  X  109 

0.1 

77 

12 

16 

0.37 

70 

7.5  X  109 

0.17 

77 

10 

25 

0.32 

64 

8X  109 

0.18 

77 

12 

10 

0.19 

18 

3.7X  109 

8 

0.43 

32 

9* 

6. 8  X  10 

*  After  antireflection  coating  with  Se. 

Table  1-1  gives  representative  values  of  responsivity,  detectivity,  and  efficiency  that  have 
been  obtained.  In  the  case  of  the  last  detector  in  the  Table,  an  antireflection  coating  of  Se  ( l/4 
wavelength  thickness  at  lOp)  was  applied  by  means  of  evaporation.  This  apparently  eliminates 


Section  I 


the  reflection  loss  over  most  of  the  8-  to  1 4 -jjl  wavelength  range,  thereby  producing  an  increase 


by  about  a  factor-of-two  in  the  quantum  efficiency. 


A.  R.  Calawa  T.C.  Harman 
I.  Melngailis  P.  Youtz 


B.  PHOTOCONDUCTIVITY  IN  SINGLE-CRYSTAL  Pb^Sn^e 

Photoconductivity  at  wavelengths  uptol5[ji  has  been  observed  at  77  °K  in  Bridgman-grown  and 
subsequently  annealed  crystals  of  Pb^_xSn  Te.  The  annealing  parameters  and  the  resulting  car¬ 
rier  concentrations  and  mobilities  of  the  crystals  used  for  photoconductors  are  shown  in  Table  1-2. 
The  two  compositions,  x  =  0.17  and  0.20,  are  of  particular  interest  for  detectors  because  for  the 
0.17  alloy  the  photoconductivity  at  77°K  peaks  very  close  to  10.6  (jl,  whereas  for  the  0.20  alloy 
the  peak  is  at  14  p,  i.e.,  detectors  of  this  composition  cover  the  8-  to  14 -p  atmospheric  window. 


TABLE  1-2 


ISOTHERMAL  ANNEALING  OF  BRIDGMAN-GROWN  Pb,  Sn  Te 

1  -x  x 


Sample 

Number 

X 

Sample 

Temperature 

(°C) 

Annealing 

Time 

(days) 

Sample 

Thickness 

(mm) 

Hall  Carrier 
Concentration  at  77°K 

(cm  3) 

Hall  Carrier 
Mobility  at  77°K 

(cm^/V  sec) 

66168C 

0.17 

500 

22 

0.29 

2.3X  1015P 

3. 1  X  104 

450 

38 

661613E2 

0.17 

650 

4 

0.27 

5.8  X  1 01 5 P 

2.6  X  104 

450 

50 

6741 

0.20 

600 

8 

0.25 

7.5  X  1015N 

3.2  X  104 

450 

56 

673R 

0.17 

600 

5 

0.06 

2.6  X  1016N 

2. 9  X 104 

450 

14 

Contacts  to  the  photoconductors  were  made  by  gold-plating  and  then  cold -bonding  leads  to 
the  gold  layer  using  indium.  The  samples  were  then  etched  both  to  reduce  their  thickness  and 

to  lower  the  surface  recombination  velocity.  The  etch  con¬ 
sisted  of  10  parts  potassium  hydroxide,  1  part  hydrogen 
peroxide,  and  10  parts  ethylene  glycol. 

In  order  to  insure  uniform  cooling  and  to  avoid  the 
danger  of  thermal  stresses,  the  samples  were  immersed 
directly  in  the  liquid  coolant  in  preference  to  mounting  them 
on  a  heatsink.  To  insure  that  photoconductivity  rather  than 
slow  thermal  effects  (bolometer  effect  or  thermoelectric 
effect)  are  observed,  the  response  speed  was  measured  by 
using  a  GaAs  diode  laser.  Figure  1-1  shows  the  response 
at  77 °K.  The  15-nsec  response  time  corresponds  to  the 
effective  carrier  lifetime.  In  addition  to  photoconductivity. 


^ —  TIME 

Fig.  1-1.  Response  of  Pbo.83^nO.  17^e 
photoconductor  at  77°K  to  GaAs  diode 
laser  pulse. 
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a  photovoltaic  effect  is  observed  at  the  gold  contacts,  presumably  due  to  the  presence  of  a 
Schottky  barrier.  This  can  be  eliminated  by  focusing  the  light  spot  on  regions  removed  from 
the  contacts. 


Fig.  1-2.  Responsivity  spectra  of  Pb]_xSnxTe 
photoconductors  at  77° K  (samples  94A  and 
97A). 


WAVELENGTH  (p) 


Fig.  1-3,  Responsivity  spectrum  of  Pbo.  83SnO.  1  7Te 
photoconductor  at  4.  2°K  (sample  94C). 


Figure  1-2  shows  the  response  spectra  of  two  photoconductors  at  77°K  —  one  for  a  0.17  and 
one  for  a  0.20  alloy.  The  response  peaks  are  at  10  and  14  p,  respectively.  Near  the  bandgap, 
there  is  a  small  peak  as  the  radiation  begins  to  penetrate  into  the  bulk  where  the  lifetime  of 
carriers  is  longer  than  near  the  surface.  The  fact  that  the  peak  is  small  indicates  that  the  sur¬ 
face  recombination  velocity  is  not  very  high.  This  means  that  the  effective  lifetime  obtained 
from  the  pulse  measurements  provides  an  approximate  measure  of  the  bulk  lifetime.  Figure  1-3 
shows  the  spectrum  of  a  PbQ  83SnQ  1?Te  photoconductor  at  4.2°K. 

Measurements  of  noise  indicate  that  for  some  types  of  contacts  there  is  a  considerable 
current  noise  with  a  1/f  type  of  frequency  dependence.  So  far,  we  have  found  that  gold-plated 
contacts  produce  the  least  noise.  For  the  photoconductors  used  at  bias  currents  up  to  the  point 
where  there  was  evidence  of  heating,  the  measured  noise  was  within  5  0  percent  of  the  Johnson 
noise  calculated  for  the  resistance  of  the  photoconductor  at  liquid  nitrogen  temperature.  The 
generation-recombination  noise  was  estimated  to  be  about  one  order  of  magnitude  smaller  than 
the  Johnson  noise  for  the  samples  used.  Hence,  Johnson  noise  appears  to  be  the  basic  noise 
limitation  at  present. 

The  properties  of  a  number  of  photoconductors  are  shown  in  Table  1-3.  At  the  900-Hz 

frequency,  the  noise  in  the  first  sample  was  about  5  0  percent  higher  than  the  Johnson  noise  for 

the  42-ohm  resistance  at  77°K.  At  3000Hz,  the  measured  noise  exceeded  the  Johnson  noise  by 

only  25  percent.  The  detectivity  is  given  for  11  p  and  900  Hz.  The  peak  D?  at  high  frequencies, 

8  l/2  ^ 

which  is  limited  by  thermal  noise,  is  about  4.5  X  10  cm/W  sec  '  . 

At  4. 2° K,  the  sensitivity  is  greatly  improved  because  the  lifetime  at  4.2°  is  much  longer 
than  at  77 °K.  Lifetime  measurements  obtained  by  observing  directly  the  decay  of  the  photo¬ 
conductivity  current  following  a  GaAs  diode  laser  pulse  are  shown  in  Table  1-4  for  a  number 
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TABLE  1-3 

PROPERTIES  OF  Pb.  Sn  Te  PHOTOCONDUCTORS 

1  -x  X 

Sample 

No. 

Crystal 

No. 

Dimensions 

Length 

Width 

Thickness 

(mm) 

T 

(°K) 

R 

(ohms) 

1 

o 

(mA) 

ft  (X  ) 

V  p 

(V/W) 

X 

P 

(h) 

D*(X  ,  900,  1) 

P 

(cm/W  sec  ^^) 

94A 

66-16 

4.0 

2.0 

0.05 

77 

42 

10 

0.7 

10 

3X  108 

94B 

66-16 

4.0 

2.0 

0.03 

4.2 

52 

4 

80 

14 

1.7X  1010 

97A 

67-4 

1.5 

1.0 

0.05 

77 

13 

30 

0.6 

14 

1  X  io8 

TABLE  1-4 

LIFETIME  MEASUREMENTS  IN  Pb,  Sn  Te 

1  -x  X 

PHOTOCONDUCTORS 

Teff 

(sec) 

T  (°K) 

Sample  94B 

Sample  97A 

300 

4  X  10"8 

— 

198 

3X  10"8 

~2  X  10_8 

77 

1.5  X 10'8 

1.5X  10"8 

4.2 

1.2  X  10"6 

0.  7  X  10-6 
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of  temperatures.  For  the  crystals  for  which  both  responsivity  and  lifetime  measurements  were 
available,  the  lifetime  calculated  from  the  responsivity  was  in  good  agreement  with  the  meas¬ 
ured  lifetime. 

The  direct  radiative  recombination  lifetime  at  77 °K  is  estimated  to  be  of  the  order  of 


10 


-5 


sec  in  Pbj_xSnxTe  with  an  energy  gap  of  0.1  eV  and  with  a  carrier  concentration  of  the 
15  -3 

order  of  10  cm  ;  hence,  additional  recombination  mechanisms  which  reduce  the  lifetime 
must  be  predominant  in  these  crystals.  ^  Melngailis 

T.  C.  Harman 


C.  LONG  WAVELENGTH  INFRARED  Pb^Sn^Te  DIODE  LASERS 

The  long  wavelength  limit  for  semiconductor  lasers  has  been  extended  to  28  p  using 

Pbi  ^Sn^Te  diode  lasers  with  x  up  to  0.27.  The  p-n  junctions  were  diffused  into  the  vapor- 

grown  single  crystals  in  a  special  fused  silica  ampoule  which  remained  sealed  throughout  the 

entire  growth  and  diffusion-annealing  processes.  Pb^  xSnxTe  diode  lasers  have  previously 

been  fabricated  by  a  different  method  for  x  4  0.17  (see  Ref.  1).  The  previous  maximum  wave 

2 

length  for  semiconductor  lasers  was  22.5  p,  observed  in  a  pressure-tuned  PbSe  laser,  while 
the  long  wavelength  limit  for  unstressed  semiconductor  lasers  was  18.8  p,  obtained  with 
Pb,  Sn  Se  (Ref.  3).  Threshold  current  densities  in  the  present  lasers  were  observed  to  be 


dependent  on  diode  surface  conditions  and  could 
inal  values  by  etching. 

A  schematic  diagram  of  the  special  growth 
and  diffusion-annealing  ampoule  is  shown  in 
Fig.I-4(a).  The  fused  silica  ampoule  consists 
of  a  tubular  section  about  1 .8  cm  i.  d.  and  14  cm 
long,  tapered  to  a  point  at  the  bottom,  and  with 
three  indentations  about  4  cm  from  the  tip.  The 

_7 

evacuated  (10  torr)  ampoule  is  suspended  in 
a  vertical  tubular  furnace,  which  has  the  tem¬ 
perature  profile  depicted  in  Fig.I-4-(b).  The 
temperature  difference  between  the  source  ma¬ 
terial  and  vapor-grown  crystals  was  about  1°C. 
The  vapor  source,  which  is  a  two-phased  coarse 
powder  obtained  from  a  quenched  ingot  of  com¬ 
position  (Pb-Sn)  q  5|TeQ  ^ ,  is  contained  in  the 
inner  quartz  tube  which  has  a  1.0  cm  i.d.  and 
is  6.5  cm  long.  The  crystals  usually  grow  on 
the  walls  of  the  ampoule  in  the  conical  tip  re¬ 
gion.  After  a  growing  time  of  approximately 

20  hours,  the  tube  is  air-cooled.  Crystals  with 
2 

(100)  surfaces  1mm  or  larger  in  area  and  of 
composition  up  to  x  =  0.27  have  been  obtained 
by  this  method.  The  rapidly  cooled,  as-grown 
crystals  have  bulk  interior  carrier  concentra¬ 
tions  which  are  given  approximately  by  the 


be  reduced  to  50  percent  or  less  of  their  orig- 


(a)  (b) 

Fig.  1-4.  (a)  Schematic  diagram  of  special  vapor 
growth  and  diffusion-annealing  ampoule  (there 
are  three  indentations  in  outer  quartz  tube  which 
are  spaced  120°  apart),  (b)  Section  of  furnace 
temperature  profile  (overall  length  of  vertically 
positioned  furnace  is  60  cm). 
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TABLE  1-5 

ANNEALING  PARAMETERS  AND  OPERATING  CHARACTERISTICS 

OF  Pb  Sn  Te  DIODE  LASERS 
1  -x  x 


X 

Anneal  ing 
Time 
(days) 

Anneal  ing 
Temperature 
(°C) 

Emission  Wavelength  (p) 

Minimum  Threshold  Current 

2 

Density  (A/cm  ) 

12°K 

77°K 

12°K 

77°K 

0.15 

None 

None 

11.7 

9.5 

250 

3,000 

0.17 

2 

450 

13.0 

9.9 

55 

10,000 

0.19 

3 

450 

14.5 

11.0 

130 

9,000 

0.20 

3 

450 

15.1 

11.2 

275 

12,500 

0.21 

2 

450 

16.8 

12.0 

175 

7,000 

0.22 

7 

450 

17.2 

- 

175 

>30,000 

0.24 

14 

400 

20.0 

- 

230 

>30,000 

0.25t 

3 

450 

22.4 

14.9 

75 

5,000 

0.27 

See  text 

See  text 

28.0 

- 

125 

>30,000 

f  The  quenched  ingot  used  for  the  vapor  growth  source  was  doped  with  4X  10 
Bi  atoms/cm3. 
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points  on  the  metal-saturated  solidus  lines  in  Ref .  4  corresponding  to  the  growth  temperature 
of  about  825  °C.  N-type  layers  of  thickness  ^25  \±  were  diffused  by  isothermally  annealing  the 
crystals  in  a  horizontal  furnace  in  the  unopened  growth  ampoule  and  thus  in  the  presence  of  the 

4 

remaining  source  ingot.  This  method  of  preparing  junctions  utilizes  the  property  that,  at  least 
for  values  of  x  up  to  x  =  0.27,  metal-saturated  crystals  change  from  p-  to  n-type  conductivity 
as  the  temperature  is  decreased.  From  isothermal  metal-saturation  experiments,  type  con¬ 
version  temperatures  of  580°,  525°,  and  425°C  are  estimated  for  x  =  0.17,  0.20,  and  0.27, 
respectively.  For  x  <  0.27,  annealing  parameters  for  junction  formation  are  given  in  Table  1-5. 
For  x  =  0.27,  the  annealing  schedule  for  preparing  junctions  was  4  hours  at  700°C,  6  days  at 
650°C,  7  days  at  450°C,  and  21  days  at  400°C,  giving  an  n-region  approximately  30  p  deep. 

Annealing  studies  indicate  the  surface  electron  concentration  of  annealed  lindoped  Pbn  77Snn  97Te 
16-3  v .  ( ~>  v .  c.  ( 

to  be  -^10  cm  .  From  the  annealing  schedules  followed  in  the  present  experiments,  we  expect 

that  the  hole  concentration  in  the  p-type  substrate  is  less  than  the  surface  n-type  carrier  density 
for  a  distance  beyond  the  n-p  junction  essentially  equal  to  the  junction  depth.  Hence,  p-n  junc¬ 
tion  phenomena  in  these  diodes  may  occur  in  a  region  of  relatively  low  carrier  concentration. 

After  annealing,  laser  structures  in  the  shape  of  rectangular  parallelpipeds  were  obtained 
by  cleaving  perpendicular  and  parallel  to  natural  (100)  faces.  To  reduce  planarity  and  reflec¬ 
tivity  of  the  longer  pair  of  sides  perpendicular  to  the  p-n  junction  plane,  these  were  subjected 
to  a  chemical  etch  consisting  of  equal  volumes  of  a  saturated  water  solution  of  Na^S^Og  and  a 
15  wt.  percent  water  solution  of  NaOH  (Ref.  5).  This  leaves  the  smaller  end  faces  to  define  the 
Fabry-Perot  cavity.  Linear  dimensions  of  the  lasers  were  typically  a  few  hundred  microns. 
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Electrical  contacts  to  the  entire  top  and  bottom  surfaces  were  made  with  electroplated  Au  fol¬ 
lowed  by  electroplated  In.  The  laser  structures  were  mounted  in  low-inductance  packages  and 
leads  were  attached  by  room-temperature  bonding  methods. 

It  has  been  found  consistently  that  a  final  fabrication  step  consisting  of  the  removal  by  an 
electrolytic  etch^  of  approximately  lOp  of  material  from  the  four  exposed  laser  sides  reduces 
the  threshold  currents  by  at  least  50  percent.  This  may  be  due  to  the  removal  of  a  damaged 

layer  at  the  cleaved  faces  which  exhibits  a  high  rate  of  nonradiative  recombination.  Etching 

7 

studies  have  shown  that  such  damage  exists  to  a  depth  of  at  least  several  microns.  This  layer 
would  decrease  the  effective  quantum  efficiency  for  recombination  of  injected  carriers  and,  by 
locally  reducing  the  density  of  injected  carriers  at  a  given  current  level,  would  produce  a  region 
of  electromagnetic  loss  near  each  end  face  of  the  laser  cavity.  The  width  of  the  lossy  region 
would  equal  approximately  the  width  of  the  damaged  layer  plus  a  diffusion  length  for  bulk  minor¬ 
ity  carriers.  Measurements  of  photoconductivity  in  Pb,  Sn  Te  at  4.2°K  indicate  that  some 

1  —x  x  g 

crystals  of  low  carrier  concentration  have  carrier  diffusion  lengths  greater  than  100p,  so  that 
even  a  narrow  high  recombination  rate  region  could  be  quite  important.  The  volt-ampere  char¬ 
acteristics  of  the  diode  lasers  exhibit  a  marked  increase  in  zero-bias  impedance  after  the  elec¬ 
trolytic  etch,  which  is  evidence  of  the  removal  of  a  region  with  a  high  recombination  rate.  In 

-3  2 

diodes  with  areas  of  2  X  10  cm  ,  etching  has  resulted  in  an  increase  of  zero-bias  impedance 
4  3 

to  as  high  as  2  X  10  ohms  for  x  =  0.20  at  12°K,  and  2X10  ohms  for  x  =  0.15  at  77°K,  from 

typical  pre-etch  values  of  5  ohms  at  both  temperatures. 

Several  lasers  were  fabricated  at  each  of  a  number  of  compositions,  and  were  evaluated  at 

77°  and  12°K.  The  experimental  arrangement  for  measuring  spectra  has  been  described  else- 
1 

where.  Emission  wavelengths  and  minimum  threshold  current  densities  at  each  composition 
are  listed  in  Table  1-5.  In  every  case,  laser  action  was  evidenced  by  a  spectrum  with  a  well- 
defined  mode  structure  exhibiting  the  mode  spacing  expected  for  the  Fabry-Perot  cavity.  The 
emission  at  28  p  by  the  PbQ  .^Sn^  ^Te  laser  sets  a  new  long  wavelength  limit  for  semiconductor 
lasers.  No  consistent  variation  of  threshold  current  with  the  Pb:Sn  ratio  is  apparent  at  either 
12°  or  77 °K.  Power  measurements  were  performed  on  a  Pb^  g^Sn^  ^Te  ^aser  using  a  cali¬ 
brated  thermopile  positioned  to  intercept  nearly  all  the  radiation  emitted  from  one  laser  end 
face.  With  pulse  bias,  a  peak  power  output  of  0.1  W  at  15  A,  corresponding  to  an  external  quan¬ 
tum  efficiency  77  of  0.08,  was  observed  at  12 °K,  and  a  peak  power  of  0.03  W  at  25  A,  giving 

x  -3 

V  x  =  0.01,  was  measured  at  77°K.  With  CW  operation  at  12°K,  5  x  10  W  at  2  A,  correspond¬ 
ing  to  77  =  0.03,  was  observed. 

The  observation  of  laser  action  up  to  x  =  0.27  indicates  that  the  bandgap  in  the  Pb^  Sn  Te 
alloy  system  is  probably  direct  to  at  least  this  composition.  The  variation  of  the  bandgap  with 
composition  as  determined  by  the  wavelength  of  the  laser  emission  is  consistent  with  the  variar 


tion  proposed  previously. 


J.  F.  Butler 
T.  C .  Harman 


D.  FAR-INFRARED  PHOTOCONDUCTIVITY  IN  HIGH-PURITY  EPITAXIAL  GaAs 

With  the  development  of  vapor  and  liquid  phase  epitaxial  techniques,  it  has  become  possible 
to  prepare  high-purity  layers  of  GaAs.  The  electrical  behavior  of  this  material  at  temperatures 
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between  300°  and  4°K  is  dominated  by  shallow  donor  impurities  with  ionization  energies  of  about 

1 0 

0.005  eV  in  the  purest  material.  We  report  here  the  observation  of  extrinsic  photoconductivity 
involving  these  shallow  impurity  levels.  The  high  responsivities  and  short  response  times  ob¬ 
served  indicate  that  high-purity  GaAs  has  considerable  potential  as  a  sensitive  and  fast  far- 
infrared  detector. 

The  results  reported  here  were  obtained  with  vapor  phase  epitaxial  GaAs  which  was  grown 
on  <11 1>  oriented  Cr-doped  semi-insulating  GaAs  substrates  using  an  AsCl^-Ga-H^  flow  sys¬ 
tem.11  The  epitaxial  layer  was  40  p  thick,  and  the  results  of  Hall  effect  measurements  gave  the 

net  donor  concentration,  Nn  —  NA,  as  6.8  x  1014cm  *  and  a  Hall  mobility  at  77 °K  of  7  5,000  cm^/ 

u  A  .  11-3 

V  sec.  The  net  carrier  concentration  at  4.2 °K  was  1.4  x  10  cm  ,  and  the  corresponding  mo¬ 
bility  was  2  500  cm2/V  sec.  The  behavior  of  the  resistivity  and  Hall  constant  below  about  8°K 
indicated  effects  of  some  impurity  banding.  Ohmic  contacts  were  formed  by  alloying  pure  Sn 
with  the  GaAs  at  300°C  in  a  reducing  atmosphere  on  two  opposite  edges  of  the  0.2  X  0.2  cm  pho- 
toconductive  samples.  The  4.2°  dark  resistance  of  these  samples  was  greater  than  3  megohms. 
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Fig.  1-5.  Voltage-current  characteristic  and  variation  of  photoconductive  responsivity 
at  337  p  and  noise  voltage  (f  =  300  Hz,  Af  =  1  Hz)  with  sample  current. 


Figure  1-5  shows  the  V-I  characteristic  of  a  photoconductive  sample  at  4.2°K.  Similar  be¬ 
havior  was  observed  in  the  other  samples,  and  the  characteristics  are  identical  in  the  reverse 
direction,  as  expected  for  ohmic  contacts.  This  behavior  is  ohmic  for  low  currents  and  voltages, 
but  as  the  electric  field  is  increased  above  about  1  V/cm,  the  V-I  characteristic  becomes  dis¬ 
tinctly  nonlinear,  with  the  current  increasing  rapidly  with  increasing  field.  At  a  critical  field 
strength  of  about  5  v/cm,  the  electric  field  decreases  to  a  lower  sustaining  value  which  remains 
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constant  with  increasing  current.  Such  behavior,  as  first  observed  in  compensated  Ge  (see 
Ref.  12),  is  due  to  impact  ionization  of  electrons  from  the  compensated  shallow  impurity  levels 
or  from  an  impurity  band;  it  has  previously  been  reported  in  GaAs  (Ref.  13). 

The  photoconductivity  was  examined  at  wavelengths  of  195  and  337  p  using  the  emission  from 

a  DCN  and  HCN  laser,  respectively,  and  at  902  p  using  the  output  of  a  carcinotron.  The  power 

level  of  the  sources  at  902  and  337  p  was  130  pW,  and  at  195  p  was  90  pW.  The  sample  was 

14 

mounted  in  a  detector  dewar  similar  to  that  used  for  submillimeter  InSb  detectors,  in  which 
the  sample  was  suspended  at  the  open  end  of  a  silver-plated  stainless-steel  light  pipe  about 
1  meter  long  and  was  immersed  in  liquid  He.  Two  black  polyethylene  filters,  one  of  which  was 
cold,  were  mounted  in  the  light  pipe  to  exclude  radiation  shorter  than  about  40  p.  A  copper  cone 
and  90°  bend  were  mounted  at  the  top  of  the  light  pipe  to  collect  the  radiation  from  the  source 
after  it  traversed  a  short  air  path.  The  photoconductive  responsivity  measured  at  337  p  is 
shown  in  Fig.  1-5  as  a  function  of  sample  current.  The  signal  increases  linearly  with  current, 
as  expected,  until  the  field  reaches  the  value  where  impact  ionization  begins  to  occur.  The 
responsivities  at  the  other  two  wavelengths  show  the  same  variation  with  current.  The  peak 
responsivities  are  5  X  10^  V/W  at  195  p,  4  x  104  V/W  at  337  p,  and  9  x  10^  V/W  at  902  p. 

The  sample  noise  was  measured  at  300  Hz  (the  chopping  frequency  used  for  the  photoconduc¬ 
tive  measurements)  with  a  low -noise,  high-impedance  preamplifier  and  a  wave  analyzer.  The 
results  of  these  measurements  are  also  shown  in  Fig.  1-5.  For  electric  fields  below  0.5  V/cm, 
the  noise  is  dominated  by  the  Johnson  noise  of  the  sample  and  the  load  resistor  (at  4.2°K)  in 
parallel  plus  a  component  due  to  the  amplifier.  However,  at  larger  electric  fields  in  the  non¬ 
linear  region  of  the  V-I  characteristic,  the  total  noise  increases  rapidly  and  becomes  dominated 
by  G-R  noise.  These  noise  measurements  can  be  combined  with  the  responsivities  to  calculate 

an  NEP.  =  V  T/R  ,  where  V  T  is  the  rms  noise  voltage,  and  R  is  the  responsivity  at  the  wave- 
AiMA  JN  -11  ^  -12 

length  A.  The  peak  values  of  NEP^  are  1.2X10  Watl95p,  1.4X10  Wat  337  p,  and 

6  x  10'11  W  at  902  (i. 

The  photoconductive  time  constant  was  determined  to  be  shorter  than  1  psec  by  using  the 

free  carrier  absorption  induced  by  impact  ionization  in  uncompensated  Ge  at  4.2 °K  to  chop  the 

incident  radiation.  A  time  constant  of  about  5  nsec  for  the  lifetime  of  the  carriers  resulting 

from  impact  ionization  of  the  shallow  impurity  levels  in  GaAs  at  low  electric  fields  was  measured 
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using  a  method  similar  to  that  described  by  Koenig. 

There  are  three  possible  mechanisms  which  must  be  considered  in  the  operation  of  a  detector 
of  this  type.  The  first  is  the  heating  of  the  sample  by  the  incident  radiation,  or  the  bolometric 
effect.  However,  since  this  effect  is  inherently  slow,  the  high  responsivities  and  short  time 
constants  observed  indicate  that  this  effect  is  not  predominant.  The  second  mechanism  is  the 
photoionization  of  the  shallow  impurity  levels  by  the  incident  radiation.  There  should  be,  a  def¬ 
inite  long-wavelength  threshold  for  this  effect  corresponding  to  the  ionization  energy  of  the  im¬ 
purity  state.  For  highest  purity  material,  this  should  occur  at  about  250 p  corresponding  to  a 
donor  level  of  ~ 0.005  eV,  but  for  material  with  higher  donor  concentrations  the  threshold  could 
occur  at  wavelengths  longer  than  400  p  corresponding  to  a  donor  ionization  energy  £0.003eV 
(see  Ref.  10).  Although  this  threshold  has  not  yet  been  located  precisely,  the  fact  that  the  photo¬ 
conductive  responsivity  at  902  p  is  much  lower  than  that  at  337  p  indicates  that  the  threshold  is 
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probably  between  these  two  wavelengths  and  that  impurity  level  photoionization  is  likely  to  be 
the  process  responsible  for  the  signal  at  the  two  shorter  wavelengths.  The  last  mechanism, 
which  is  probably  responsible  for  the  response  at  902  p,  is  the  free  carrier  hot  electron  effect 
as  observed  in  InSb  (see  Ref.  16).  The  coupling  between  the*  electrons  and  the  crystal  lattice  is 
very  weak  at  low  temperatures  in  these  materials  so  that,  even  at  very  low  electric  fields,  the 
average  energy  of  the  photo-excited  carriers  can  be  appreciably  greater  than  the  thermal  equilib¬ 
rium  value.  Since  the  mobility  is  limited  by  ionized  impurity  scattering  at  low  temperatures, 
the  mobility  of  the  hot  electrons  is  increased  and  the*  sample  resistance  decreases.  The  response 
due  to  this  mechanism  should  be  most  important  at  long  wavelengths. 

G.E.  Stillman  C.  D.  Parker 

C.M.  Wolfe  P.E.  Tannenwald 

1.  Melngailis  J.O.  Dimmock 

E.  ISOLATION  OF  JUNCTION  DEVICES  IN  GaAs  USING  PROTON 
RADIATION  DAMAGE 

We  have  found  that  proton  radiation  damage  can  be  used  to  convert  both  p-  and  n-type  GaAs 
into  high-resistivity  material,  and  have  used  this  technique  to  isolate  p-n  junctions  on  a  diffused 
GaAs  substrate,  to  isolate  islands  of  n-type  material  in  an  epitaxial  layer  deposited  on  semi- 
insulating  GaAs,  and  also  to  prevent  edge  breakdown  in  Au-GaAs  Schottky  barriers. 

The  ion  generator  used  in  these  experiments  is  a  modified  400-keV  Van  de  Graaff  generator 
designed  to  accelerate  positive  ions.  The  protons  are  generated  by  RF  excitation,  accelerated 
to  400  keV,  and  magnetically  analyzed.  The  samples  are  mounted  on  a  tantalum  sample  holder, 
with  part  of  each  sample  shielded  from  the  ion  beam  with  a  molybdenum  mask. 

The  p-n  junction  devices  were  formed  on  a  slice  of  n-type  GaAs  with  n  =  2  x  10^/cm^. 
Following  the  diffusion  of  a  Zn-doped  p-type  layer  into  the  top  surface,  an  array  of  0.003-inch 
gold  circles  on  0. 005-inch  centers  was  plated  onto  the  sample  using  0.1  5-p  pyrolytic  SiO^  as  a 
plating  mask.  These  gold  contacts  were  sufficiently  thick  (3  p)  that  the  proton  beam  could  not 

o 

penetrate  them.  However,  the  1500-A  SiO^  layer  between  the  contacts,  which  was  left  on  the 

sample  to  reduce  sputtering  effects,  was  sufficiently  thin  so  that  it  had  only  a  slight  effect  on 

the  proton  beam.  Consequently,  only  that  part  of  the  sample  between  the  contacts  was  bombarded. 

Each  sample  was  then  irradiated  at  room  temperature  with  a  proton  beam  of  current  density 

0.1  (lA/cm^,  starting  with  an  energy  of  400  keV  and  subsequently  decreasing  the  beam  energy  in 

steps  of  20keV,  finishing  at  60keV.  At  each  energy,  the  sample  was  exposed  for  20  sec,  to  give 
14  2 

a  total  dose  of  about  10  protons/cm  .  Before  irradiation,  the  current-voltage  characteristics 
between  the  contacts  were  linear  with  a  resistance  of  about  500  ohms.  Following  irradiation, 
these  same  contacts  gave  back-to-back  diode  characteristics,  with  a  leakage  current  less  than 
0.1  fjiA  and  a  breakdown  voltage  of  about  40  V,  typical  of  the  substrate.  Between  a  gold  contact 
on  the  p-layer  and  an  ohmic  contact  on  the  substrate,  normal  diode  characteristics  were  observed. 

The  islands  of  n-type  material  were  prepared  on  a  sample  which  had  a  3-|j.  thick  epitaxial 
layer  of  n-type  material  (n  ~  10^/cc)  which  was  deposited  on  a  chromium-doped  semi-insulating 
substrate.  The  sample  was  plated  with  an  array  of  gold-tin  contacts  which  were  alloyed.  Ad¬ 
ditional  gold  was  plated  to  provide  a  thickness  of  3  p.  The  sample  was  then  irradiated  at  400  keV 
for  20  sec.  Before  irradiation,  the  contact-contact  resistance  was  about  60ohms;  subsequently. 
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the  leakage  currents  were  less  than  0. 1  (jlA  at  100  V  between  these  same  contacts,  indicating 

q 

that  the  leakage  resistance  between  the  islands  was  in  excess  of  10  ohms. 

The  Schottky  barrier  diodes  were  prepared  in  a  manner  similar  to  the  p-n  junction  devices, 
except  that  no  p-type  layer  was  diffused  into  the  substrate,  and  irradiation  was  done  only  at 
400  keV.  Following  irradiation,  the  original  3-p  thick  gold  contacts  were  removed  and  slightly 
larger  semitransparent  gold  contacts  were  plated  on  the  same  areas.  Current-voltage  measure¬ 
ments  on  these  devices  showed  that  the  reverse  breakdown  voltage  following  irradiation  had 
increased  substantially  compared  with  the  original  diode,  while  the  forward  characteristic  was 
unchanged.  Visual  observation  of  the  light  emitted  from  these  diodes  in  reverse  bias  breakdown 
indicated  that  edge  breakdown  which  was  present  in  the  original  diodes  had  been  suppressed,  and 
that  the  higher  voltage  breakdown  in  the  irradiated  diodes  occurred  away  from  the  edge. 

We  have  also  done  some  experiments  in  which  the  sample  was  bombarded  before  the  gold 

contacts  were  applied,  in  order  to  determine  the  doping  profile  of  the  irradiated  regions.  These 

experiments  were  done  by  bombarding  samples  of  n-type  GaAs  with  proton  beams  of  intensity 

0.1  pA  for  20  sec  at  various  energy  levels.  After  bombardment,  each  sample  was  plated  with 

gold  contacts,  and  capacity-voltage  measurements  were  done  in  order  to  determine  the  doping 

profile.  The  initial  results  of  these  experiments  are  that  the  damaged  region  has  high  resistance 
1 3 

(doping  levels  ~10  /cc  or  less)  from  the  surface  of  the  crystal  until  the  protons  approach  the 
end  of  their  range.  Near  the  end  of  the  range,  there  is  evidence  of  a  change  in  the  character 
of  the  radiation  damage,  with  an  indication  that  shallow  donors  are  created  in  this  region.  The 
depth  of  the  high-resistance  region  is  about  4  p  at  400  keV. 

This  technique  of  creating  high-resistivity  material  is  simple  and  quite  easy  to  use.  We 
feel  that  it  should  be  a  useful  tool  for  the  fabrication  of  other  structures  in  which  thin  high- 


resistivity  layers  are  required. 


A.G.  Foyt  J.  P.  Donnelly 

W.T.  Lindley  C.M.  Wolfe 


F.  SHORT  WAVELENGTH  RESPONSE  CHARACTERISTICS 
OF  InSb-MOS  PHOTOVOLTAIC  DETECTORS 

We  have  continued  to  study  the  effects  of  short  wavelength  light  on  the  characteristics  of 
17 

InSb-MOS  detectors.  The  photovoltaic  response  to  modulated  3.9-^  radiation  from  an  InAs 
diode  emitter  has  been  measured  at  77 °K.  In  our  model  of  the  MOS  structure,  this  radiation 
acts  as  a  probe  of  the  InSb  depletion  region,  and  has  no  effect  on  the  state  of  the  structure.  For 
n-type  InSb,  the  depletion  layer,  in  turn,  gives  an  indication  of  the  charge  state  of  the  oxide. 

The  surface  of  the  oxide  is  illuminated  with  short  wavelength  radiation  from  a  monochromator, 
and  the  change  in  the  response  to  the  3.9-p  radiation  is  measured.  At  77 °K,  there  is  no  response 
to  the  InAs  radiation  so  long  as  the  constant  irradiation  has  a  wavelength  corresponding  to  photon 
energies  <  1 . 3  eV  or  >4,8eV.  In  equilibrium  for  photon  energies  >1 . 3  eV  and<4.8eV,  the  re¬ 
sponse  to  the  3.9-p  radiation  is  relatively  constant,  independent  of  the  intensity  of  the  constant 
irradiation,  but  the  time  to  reach  the  final  value  depends  linearly  on  the  intensity  of  the  irra¬ 
diating  light,  and  very  strongly  on  the  photon  energy. 

At  liquid  nitrogen  temperatures,  once  a  response  to  the  3.9-p  radiation  is  created  by  short 
wavelength  irradiation,  the  response  remains  for  at  least  several  hours,  even  if  the  irradiating 
light  is  cut  off.  In  fact,  the  response  remains  even  if  the  MOS  structure  is  warmed  to  room 


11 


Section  I 


temperature  and  then  cooled  again.  However,  we  find  that  if  the  sample  is  warmed  to  room 
temperature  and  then  irradiated  with  short  wavelength  photons  before  it  is  cooled  for  the  next 
measurement,  the  response  to  the  3.9-p  radiation  is  eliminated.  The  best  wavelength  for 
'’killing'*  the  response  seems  to  occur  at  the  peak  of  the  oxide  response  shown  in  Fig.  1-4  of 
Solid  State  Research  Report  (1968:1).  These  effects  are  undoubtedly  due  to  trapping  and  detrap¬ 
ping  of  electronic  charge  at  oxide  states  and/or  at  oxide-InSb  interface  states;  however,  the  pre¬ 
cise  model  for  this  trapping  has  not  yet  been  determined. 

W.  E.  Krag 

R.  J.  Phelan,  J r. 

J .  O.  Dimmock 
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n.  OPTICAL  TECHNIQUES  AND  DEVICES 


A.  100-WATT  C02  AMPLIFIER 

Further  measurements  have  been  taken  of  gain,  beam  size,  and  frequency  pulling  on  the 
9-meter  100-W  CO^  amplifier.  Details  of  the  device  are  given  in  the  previous  Solid  State 
Research  Report.! 

Gain  measurements  were  made  with  a  sealed-off  gas  mixture  of  CC^,  N^.  He,  and  H^  in 
the  amplifier  and,  although  not  optimized,  the  gain  was  within  30  percent  of  that  obtained  with 
a  flowing  gas  mixture.  This  probably  indicates  that  under  the  appropriate  conditions  the  am¬ 
plifier  gain  would  be  the  same  for  flowing  or  sealed-off  gas  mixtures.  Detailed  gain  measure¬ 
ments  for  each  tube  show  that  the  power  gain  was  fairly  well  approximated  by  a  solution  to  the 
gain  equation  assuming  a  constant  gain  with  radial  position  and  the  power  gain  only  a  function 
of  amplifier  length  L.  The  homogeneous  saturation  equation  is 


dP 

dz 


a 

o 


P 

P 

1  +  w  s 

O 


and,  when  solved  under  the  restrictions  above,  yields  (when  integrated)  the  gain  equation 
P. 

InG  +  (G  -  1)  =  a  L 

Ps 

where 

G  =  P/P.^,  the  gain  at  L  within  the  amplifier, 

P.^  =  the  input  power, 

P  =  saturation  power  =  W  S,  the  beam  cross  section  is  S,  and 
s  o 

- 1 

c*q  =  cm  ,  the  small  signal  gain  parameter. 

Figure  II—  1  shows  the  theoretical  gain  curves  as  a  function  of  the  ratio  of  tube  input  power  to 
saturation  power  of  the  amplifier.  The  parameters  of  these  curves  were  determined  from  the 
experimental  data  and,  over  the  region  shown,  represent  the  data  fairly  well.  The  power  output 
measured  for  2-  and  10-W  input  was  56  and  88  W  in  the  sealed-off  condition,  which  would  be  73 
and  114  W  for  a  flowing  gas  mixture.  The  small  signal  gain  for  the  20-mm  bore  tube  is  9.3,  the 
25-mm  bore  tube  is  8,  and  the  30-mm  bore  tube  is  4.2. 

The  beam  profile  as  shown  in  Fig.  II- 2  is  taken  at  the  output  of  the  9-meter  amplifier.  The 
solid  curve  is  Gaussian  and  closely  approximates  the  measured  laser  beam  radial  distribution 
of  the  TEMqq  mode  at  the  amplifier  output  when  the  amplifier  is  in  the  off  condition;  but,  with 
amplification,  this  approximation  becomes  inaccurate  with  increasing  beam  radius,  presumably 
because  saturation  effects  are  less  at  the  larger  radii,  resulting  in  a  larger  gain  and  higher  rel¬ 
ative  power  at  larger  distances  from  the  tube  center. 

t  Solid  State  Research  Report,  Lincoln  Laboratory,  M.I.T.  (1968:1),  p.  11,  DDC  668762. 
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Fig.  11-1.  CO2  amplifier  tube  gain  as  function  of  ratio  of  input  power  to 
saturation  power. 


Fig.  11-2.  Beam  profile  of  laser  and  amplifier  output  compared  with  Gaussian 
beam  shape. 
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The  output  of  the  laser  has  been  monitored  with  a  monochromator  which  will  resolve  a  fre¬ 
quency  difference  of  10  GHz  at  10  pm.  The  laser  output  within  this  resolution  limitation  is  not 
shifted  with  the  amplifier  on  or  off,  and  independent  of  a  static  or  flowing  gas  fill  condition  ex- 
isting  in  the  amplifier  tubes.  R  j  Carbone 


B.  STABLE  LASER  OSCILLATOR,  DESIGN  AND  TESTING 

The  extensive  testing  of  the  1^-meter  cavity  stable  oscillator  has  been  completed.  Since 
its  completion  over  a  month  ago,  this  laser  was  successfully  utilized  as  a  single-frequency 
driver  for  the  100-W  amplifier  described  in  Sec.  II-A  above. 

Two  potential  problems  arose  during  tests  of  this  10-  to  15-W  oscillator.  At  first,  in  spite 
of  the  100-MHz  cavity  mode  spacing,  this  laser  would  only  oscillate  in  the  00°  1  to  02° 0  vibra¬ 
tional  band  of  CO^  during  several  weeks  of  continuous  testing.  Although  systematic  variations 
of  pressure  and  excitation  current  have  been  tried  over  more  than  a  2:1  range,  the  output  would 
at  most  shift  to  some  other  strong  P-  or  R-branch  rotational  transition  within  the  00°  1  to  02 °0 
vibrational  band.  Subsequently,  this  behavior  was  attributed  to  slightly  higher  reflectivities  in 
the  9-  to  10-pm  wavelength  region  on  several  output  mirrors  used  in  these  tests.  After  changing 
to  mirrors  with  reflectivities  carefully  maximized  at  10.6  pm,  00 °1  to  10°0  vibrational  band 
operation  was  immediately  obtained. 

The  second  potentially  very  troublesome  problem  was  located  when  it  was  noticed  that  the 
output  mirror  holder's  temperature  increased  appreciably  with  some  of  the  mirrors  tried  during 
the  testing  of  the  1 meter  oscillator.  Subsequently,  it  was  determined  that  most  of  the  IRTRAN 
IV  substrates  tested  had  an  absorption  appreciably  greater  than  that  specified  by  the  manufac¬ 
turer.  Because  of  the  unexpected  absorption  losses  in  the  output  mirror  and  the  high  output 
power  of  the  l|- meter  oscillator,  the  output  may  shift  several  rotational  transitions  before 
reaching  a  final  equilibrium  temperature.  Although  this  settling  time  is  very  much  longer  than 
found  with  previous  lower  power  tubes,  essentially  single-transition,  single-frequency  opera¬ 
tion  will  ultimately  occur. 

The  absorption  problem  found  with  IRTRAN  IV  will  be  hopefully  overcome  by  the  future  use 
of  GaAs  substrates.  All  three  GaAs  substrates  tested  so  far  yielded  consistent  data  which  were 
well  below  the  absorption  losses  of  the  best  IRTRAN  IV  tested. 

Extensive  modifications  of  the  laser  processing  station  were  carried  out.  The  resulting 
versatility  in  the  choice  of  gas  fill  enabled  us  to  double  to  3  W  the  single-frequency  output  power 
of  the  older  50-cm  cavity  (20-cm  discharge  length)  stable  lasers.  Assembly  of  the  completely 
redesigned  versions  of  the  stable  lasers  was  begun,  and  testing  should  commence  within  a  month. 

C.  Freed 


C.  THERMAL  BLOOMING 

From  elementary  considerations,  one  would  expect  lens-like  effects  when  a  light  beam  de¬ 
posits  heat  in  a  fluid.  Assuming  the  refractive  index  change  is  proportional  to  the  density  change, 

1  2 

analytic  expressions  for  the  steady-state  focal  length  (negative)  and  the  time  development  of 
the  focal  length  of  an  absorbing  medium  can  be  obtained  using  a  simplified  model.  The  validity 
of  this  model  has  been  experimentally  confirmed  in  liquid  cells.^ 
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A  cell  has  been  constructed  to  attempt  the  measurement  of  these  thermal  blooming  effects 

in  gases.  Because  of  the  relatively  high  power  available,  a  CO^  gas  laser  has  been  used  as 

the  source  and  pure  CO^  gas  at  about  one  atmosphere  as  the  absorber.  When  necessary,  the 

absorption  can  be  enhanced  by  the  addition  of  a  few  torrs  of  ethylene.  Although  the  blooming 

in  a  2-meter  cell  with  5-W  input  in  a  1-cm  beam  is  quite  substantial  (f  ~-2m),  all  attempts  at 

quantitative  measurements  of  the  time  development  have  been  unsuccessful.  The  dynamic 

4 

range  and  sensitivity  of  the  temperature-sensitive  phosphor  are  inadequate  with  the  laser  power 
now  available  to  us  (the  beam  is  visible  with  appropriate  mounting  of  the  phosphor,  but  the  time 
response  under  these  conditions  is  poor).  An  attempt  was  made  to  measure  the  output  beam 
diameter  by  scanning  it  across  a  photoconductive  detector,  but  the  lack  of  cylindrical  symmetry 
due  to  convection  and  the  shape  changes  as  steady  state  is  approached  make  interpretation  of 
the  results  difficult.  The  cell  is  now  mounted  horizontally,  and  the  characteristic  time  for  con¬ 
vection  effects  is  comparable  to  that  for  thermal  diffusion,  causing  the  above  problems. 

Several  possible  changes  are  under  consideration,  such  as  using  a  wavelength  where  better 
image  detectors  are  available,  using  higher  power  at  10.6  p,  or  using  the  10.6  beam  only  for  heat 
deposition  and  using  a  low-power  visible  beam  for  measurements. 

A.  H.  M.  Ross 

D.  OPTICAL  HETERODYNE  DETECTION  AT  10.6  pm  OF  THE  BEAT  FREQUENCY 
BETWEEN  A  TUNABLE  PbQ  88SnQ  12Te  DIODE  LASER  AND  A  C02  GAS  LASER 

We  have  used  a  current-tunable  PbQ  g8^n0  j^e  diode  laser  to  achieve  optical  heterodyning 
with  a  CO^  gas  laser.  The  beat  frequency  was  continuously  controllable  from  under  50  to 
1300  MHz,  the  high-frequency  limit  of  our  Ge:Cu  detector.  These  experiments  have  yielded  the 
first  direct  measurement  of  the  linewidth  of  the  emission  from  a  semiconductor  p-n  junction 
laser,  displayed  on  a  spectrum  analyzer  as  a  frequency  spread  of  the  heterodyne  signal.  Pre¬ 
liminary  indications  are  that  the  linewidth  is  less  than  100  kHz,  with  power  supply  ripple  in  the 
diode  laser  circuit  presently  precluding  a  more  accurate  determination.  (The  linewidth  of  the 
CO^  laser  emission  is  less  than  a  few  hundred  hertz. 

The  diode  used  for  the  heterodyne  experiments  was  cleaved  into  a  rectangular  parallelepiped 
1.3  mm  long,  0.5  mm  wide,  and  0.2  mm  thick  from  vapor-grown  Pb^  xSnxTe  material  tailored  to 
produce  spontaneous  emission  at  10.6  pm,  for  which  x  =  0.12  Ohmic  contacts  were  applied  and 
the  diode  was  attached  to  the  cold  "finger”  of  a  liquid  helium  dewar.  A  single-pass  grating 
spectrometer  was  used  to  measure  the  wavelength  of  the  laser  emission  as  a  function  of  diode 
current;  curves  similar  to  those  of  magnetic-field  tuning^  were  obtained. 

_  i 

The  diode  laser  line  is  tunable  through  six  modes  from  943.6  to  946.7  cm  for  direct  cur¬ 
rents  from  675  to  1800  mA.  (At  1000  mA,  the  output  power  is  approximately  lOOpW. )  Within 
each  mode,  continuous  control  is  achieved  over  approximately  0.2  cm  (6  GHz).  The  overall 
tuning  range  encompasses  two  adjacent  CO^  rotational  transitions:  P^g  at  945.93  cm  ,  and 
P^q  at  944.15  cm  Only  the  P^q  transition  is  close  enough  to  one  of  the  diode  laser  modes  to 

produce  a  beat  frequency  within  the  limits  of  our  detector.  The  C09  laser,  which  was  especially 

5  z 

constructed  for  highly  stable,  single-frequency  operation,  was  tuned  to  the  P  ^  transition;  the 
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PbQ  88^n0  i2Te  diode  current  was  adjusted  to  produce  laser  emission  in  the  proper  mode.  Both 
beams  were  directed  collinearly  onto  the  Ge:Cu  detector,  and  under  these  conditions  the  hetero¬ 
dyne  IF  detector  current  is  given  by 

iIF  =  (2qGjj  /hv )  (P1P2)1/2  sin[27rfIFt]  (1) 

where  the  beat  frequency  fjp  =  c(944.15  —  and  are  the  absorbed  powers,  G  is  the 

detector  gain,  and  77  is  its  efficiency.  Polarization  and  nonalignment  terms  are  neglected. 

Figure  II- 3(a)  is  an  oscillograph  of  the  100-MHz  beat  signal  produced  between  the  P^Q  laser 
line  and  the  diode  laser  line  emitted  at  1135  mA.  The  sinusoidal  variation,  in  accordance  with 
Eq.  (1),  is  a  clear  indication  of  heterodyne  operation.  The  dispersion  of  a  1250-MHz  beat  fre¬ 
quency  (diode  current  of  1045  mA),  occurring  just  below  the  detector  cutoff  frequency,  is  shown 
in  Fig.  II- 3(b ).  It  is  evident  from  the  width  of  the  spectrum  that  the  diode  laser  linewidth  is 
100  kHz  or  less.  The  1-sec  exposure  of  Fig.  II- 3( c )  represents  100  sweeps  of  the  spectrum  an¬ 
alyzer  and  illustrates  the  fluctuations  in  the  beat  frequency  caused  by  the  100-jjlA  rms  ripple  in 
the  diode  laser  current.  With  a  better  power  supply  in  conjunction  with  storage  batteries,  we 
hope  to  reduce  the  ripple  in  the  diode  current. 

Figure  II- 4  shows  the  range  over  which  the  beat  frequency  was  continuously  tunable,  and 
the  dependence  of  the  beat  frequency  on  diode  current.  The  tuning  rate  of  12.8  MHz/mA  agrees 
well  with  that  predicted  from  the  grating  spectrometer  data. 

E.  D.  Hinkley 
T.  C.  Harman 
C.  Freed 
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Fig.  11-3.  Beat  frequency  signal  between  PbQ.  QgSng.  I2^e  diode  laser  line  and  P20  line 
of  C02  laser,  (a)  f)F  =  100 MHz,  ld  =  1135mA;  (b)  fJF  =  1250MHz,  ld  =  1045mA,  sin¬ 
gle  scan;  (c)  multiple  (100)  scan,  illustrating  effect  of  power  supply  ripple  of  100-pA 
rms  in  diode  circuit. 


Fig.  11-4.  Tuning  range  of  beat  frequency  and 
dependence  on  diode  current.  Frequency  values 
were  measured  with  spectrum  analyzer. 
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III.  MATERIALS  RESEARCH 


A.  BOUNDARY  LAYER  ANALYSIS  OF  FORCED  CONVECTION  CRYSTAL  GROWTH 

We  have  found  that  forced  convection  can  be  used  to  increase  the  rate  of  crystal  growth  from 

1-3 

the  vapor  phase  and  to  control  the  morphology  of  the  growing  interface.  The  relationship  be¬ 
tween  growth  rate  and  carrier  gas  velocity  has  now  been  derived  theoretically  by  analyzing  the 
properties  of  the  diffusion  boundary  layer  at  the  interface.  In  addition,  the  condition  required 
for  interface  stability  has  been  formulated  in  terms  of  the  temperatures  at  the  interface,  at  the 
edge  of  the  boundary  layer,  and  at  the  solid  source. 


Fig.  Ill-l.  Schematic  drawing  of  apparatus  for  measuring  crystal  growth  velocity 
as  function  of  carrier  gas  velocity  and  temperatures  T.,  T^,  and  T^. 


The  apparatus  used  to  study  forced  convection  crystal  growth  is  shown  in  Fig.  Ill—  1 .  An 
inert  carrier  gas,  passing  through  a  packed  bed  of  polycrystalline  material,  is  saturated  with 
vapor  at  the  source  temperature  Tg.  The  concentration  Cg  of  the  vapor  is  related  to  the  satu¬ 
ration  vapor  pressure  pg  by  the  ideal  gas  law  C  =  p/RT.  The  saturated  gas  is  then  heated  to  a 
higher  temperature  T^.  Finally,  it  passes  over  the  growing  crystal  interface  at  a  lower  temper¬ 
ature  T.,  where  condensation  reduces  the  vapor  pressure  to  p.  (the  saturation  pressure  at  T.), 
and  correspondingly  reduces  the  concentration  to  C..  If  it  is  assumed  that  the  temperature  and 
concentration  changes  occur  in  a  boundary  layer  of  thickness  6  adjacent  to  the  interface,  forced 
convection  may  be  regarded  as  placing  a  virtual  vapor  source  (at  temperature  T  )  and  a  virtual 
heat  source  (at  temperature  T^)  at  a  distance  6  from  the  interface. 

If  diffusion  through  the  boundary  layer  is  the  rate-limiting  step  in  crystal  growth,  the  growth 
rate  given  by  the  material  balance  at  the  interface  is 

MD(Cg-Ci) 
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where  M  is  the  molecular  weight  of  the  condensing  vapor,  D  is  its  diffusion  coefficient  in  the 
carrier  gas,  and  p  is  the  density  of  the  crystal.  For  forced  convection,  the  thickness  of  the 
boundary  layer  is 

.1/2 


6  =  B(D/V) 


(2) 


where  V  is  the  carrier  gas  velocity,  and  B  is  a  geometric  constant.  Substituting  Eq.  (2)  into 
Eq.  (1 )  gives 


M(CS  -  C.)  (DV) 

pB 


1/2 


(3) 


Thus,  the  crystal  growth  velocity  is  proportional  to  the  first  power  of  the  available  vapor  con¬ 
centration  (C  —  C^),  and  to  the  square  root  of  the  carrier  gas  velocity. 

Data  previously  reported*  for  the  growth 
velocity  of  iodine  crystals  as  a  function  of  the 
temperature  gradient  at  the  interface  can  be 
used  to  test  Eq.  (3).  These  data,  which  were 
obtained  with  the  apparatus  of  Fig.  Ill- 1  for  two 
different  helium  carrier  gas  velocities,  have 
been  replotted  in  Fig.  Ill- 2  as  a  function  of  the 
available  vapor  concentration.  In  agreement 
with  Eq.  (3),  it  is  seen  that  v  is  proportional 
to  (C  —  C.)  for  a  given  value  of  V.  The  straight 
lines  shown  in  Fig.  Ill— 2  were  calculated  from 
Eq.  (3)  by  using  the  measured  values  of  (C  —  C.) 
and  V,  the  published  values  of  M  and  p,  and 
D  =  0.5  cm  /sec,  the  self-diffusion  coefficient 
for  helium  gas  at  one  atmosphere.  The  only 

adjustable  parameter  was  the  constant  B.  The 

1  /2 

same  value  of  this  constant,  0.71  cm  '  ,  was 
used  to  calculate  the  lines  for  both  values  of  V. 
The  satisfactory  agreement  between  the  lines 
and  the  experimental  points  shows  that  the  data 
are  consistent  with  the  v  cc  v*/2  relationship  of 
Eq.  (3). 

The  experimental  results  previously  re- 
2 

ported  for  the  growth  rate  of  camphor  crystals 
also  confirm  the  proportionality  between  v  and 


WO 

Fig.  MI-2.  Growth  velocity  of  iodine  crystals  as 
function  of  available  concentration  (C$  —  Cj),  for 
two  differenthelium  carrier  gas  velocities.  Source 
temperature  Ts  =  90°C. 


(Cs-Ci). 


The  best  fit  to  the  data  was  obtained 

1/2 


for  B  =  0. 57  cm 

Throughout  the  preceding  discussion,  it 

was  assumed  that  crystal  growth  takes  place  at 
1-3 


a  smooth  interface.  However,  we  showed  earlier*  ^  that  roughening  of  the  interface  or  even 
dendrite  formation  can  occur  in  vapor  growth  because  of  interface  instability  resulting  from 
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3 

constitutional  supercooling.  We  derived  a  criterion  for  interface  stability  in  terms  of  the 
temperature  gradient  in  the  gas  phase  at  the  interface  and  the  growth  velocity,  considered  as 
two  independent  variables.  We  shall  now  extend  the  earlier  treatment  by  taking  account  of  the 
relationship  between  these  quantities. 

It  was  shown  in  Ref.  3  that  the  growth  interface  will  be  stable  if  (dp/dx).  <  (dpQ/dx)^,  where 
p  is  the  partial  pressure  of  the  condensing  component,  pQ  is  the  saturated  vapor  pressure,  and 
x  is  the  distance  between  the  interface  and  a  point  in  the  gas  phase.  The  stability  criterion  may 
also  be  written  as  (dC/dx)^4  (dCQ/dx)^,  where  CQ  is  the  concentration  corresponding  to  pQ. 

If  the  concentration  gradient  within  the  boundary  layer  is  assumed  to  be  constant,  then  for 
forced  convection, 

C  —  C. 

(dC/dx).  =  — 2-z - i  .  (4) 

1  o 

If  it  is  assumed  that  pQ  is  proportional  to  exp  [-  AH/RT]  and  that  p.  is  equal  to  pQ  at  T.,  then 
differentiating  Cq  =  pQ/RT  gives 

(dCc/dT)i  =  t:  (It  - *)  ■  <5> 

i  i 


If  it  is  assumed  that  the  temperature  gradient  is  constant  within  the  thermal  boundary  layer, 
and  that  this  layer  has  the  same  thickness  as  the  diffusion  boundary  layer,  then 


(dT/dx).  = 


(6) 


Multiplying  Eqs.  (5)  and  (6)  gives 


(dC  /dx).  = 
o'  1 


T*  -  T.  C. 

6  i  _ l 

6  T, 


(7) 


To  obtain  the  criterion  for  interface  stability,  we  set  (dCQ/dx).  equal  to  (dC/dx)..  From 
Eqs.  (4)  and  (7), 


T  -  T.  (C  -C.)/C. 
6  i  _  s  l  ’  l 

T.  =  (AH/RT.)  -  1 


(8) 


T6  (Ti/Ts)  exPKAH/RTi)  (TS/Ti“  1)1-1 
T7  =  (AH/RT.)  -  1 


(9) 


The  growth  interface  will  be  stable  only  for  values  of  T^/T^  equal  to  or  greater  than  those  given 

by  Eq.  (9).  Values  of  T./T.  calculated  from  Eq.  (9)  are  plotted  in  Fig.  Ill- 3  as  a  function  of 

T  /T.  for  the  case  in  which  AH/RT.  =  15,  the  value  for  a  gas  obeying  Trouton' s  rule  when 

p.  =  15  torrs.  We  see  from  Fig.  Ill- 3  that  along  the  calculated  line  the  value  of  T^/T.  is  always 

greater  than  the  corresponding  value  of  Tg/T^.  Therefore,  stable  growth  can  occur  only  if  the 

carrier  gas,  after  passing  over  the  source  at  T  ,  is  heated  to  a  temperature  higher  than  T  be- 

s  s 

fore  it  reaches  the  growth  interface.  It  should  be  noted  that  the  growth  velocity  does  not  appear 
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in  Eq.  (9).  Therefore,  once  a  suitable  temperature  profile  is  established,  the  growth  velocity 
may  be  increased  indefinitely,  without  causing  interface  instability,  by  increasing  the  carrier 
gas  velocity.  T  B  Reed 

A. J.  Strauss 

B.  OPTICAL  TRANSMISSION  OF  Cd^Mg^e  ALLOYS 

The  first  investigation  of  the  pseudo-binary  CdTe-MgTe  system  was  recently  reported  by 
Itoh.  Since  CdTe  has  the  zincblende  structure  and  MgTe  has  the  wurtzite  structure,  a  complete 
series  of  solid  solutions  cannot  be  formed.  However,  Itoh  found  that  there  is  extensive  mutual 
solid  solubility.  Alloys  with  zincblende  structure  are  formed  for  compositions  up  to  about 
60  mole-percent  MgTe,  and  wurtzite-structure  alloys  are  stable  for  compositions  above  about 
7  5  mole-percent  MgTe.  Itoh  found  that  the  energy  gap  of  the  zincblende-structure  alloys  in- 

4 

creases  monotonically  with  increasing  MgTe  content. 

In  the  present  study,  optical  transmission  measurements  have  been  made  on  alloy  crystals 

containing  up  to  46  mole-percent  MgTe.  The  alloys  were  prepared  from  the  elements  in  an  open- 
5 

tube  reactor,  where  synthesis  can  be  carried  out  quite  rapidly  without  the  danger  of  explosion 
encountered  when  sizable  quantities  of  the  elements  are  reacted  in  a  sealed  quartz  ampoule.  To 
obtain  single  crystals,  ingots  of  the  alloys  were  prepared  by  the  Bridgman  method.  The  fine¬ 
grained  material  prepared  in  the  reactor  was  placed  in  a  24-mm  i.d.  fused  silica  ampoule,  tapered 
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to  a  point  at  the  lower  end,  which  had  been  coated  with  pyrolytic  graphite.  After  the  ampoule  was 

evacuated  and  sealed,  it  was  heated  in  a  vertical  resistance  furnace  to  melt  the  alloy,  and  then 

lowered  through  the  furnace  to  produce  directional  crystallization  of  the  melt.  Thermal  analysis 

experiments  on  CdTe-rich  melts  have  shown  that  the  liquidus  temperature  initially  increases 

linearly  with  increasing  MgTe  content.  For  35  mole-percent  MgTe,  the  liquidus  temperature 

4 

is  1150°C.  The  value  of  1114°C  reported  by  Itoh  for  this  composition  is  much  too  low. 

Single-crystal  samples  between  1  and  3  mm  thick  were  cut  from  the  Bridgman-grown  ingots 
with  a  wire  saw  and  then  mechanically  polished.  After  the  initial  transmission  measurements, 
a  number  of  samples  were  annealed  several  days  in  saturated  Cd  vapor  at  900 °C  and  then 
remeasured. 


Fig.  IH-4.  Optical  transmission  at  room  temperature  vs  wavelength  for  Cd^^Mg^Te  alloys. 

Optical  transmission  measurements  in  the  vicinity  of  the  fundamental  absorption  edge  were 
made  at  room  temperature  with  a  single-beam,  prism  spectrometer.  Figure  II I— 4  shows  the 
transmission  curves  obtained  for  CdTe  and  three  representative  alloy  samples.  These  results 
illustrate  the  monotonic  decrease  in  the  absorption  edge  wavelength  with  increasing  MgTe  content, 
from  about  0.86  pm  in  the  infrared  for  CdTe  to  about  0.56  pm  in  the  visible  for  x  =  0.46.  They 
also  illustrate  the  marked  increase  in  transmission  at  wavelengths  beyond  the  absorption  edge 
which  generally  resulted  from  annealing  in  saturated  Cd  vapor.  (The  as-grown  data  for  the 
sample  with  x  =  0.24  have  been  omitted  from  Fig.  II I- 4 ,  since  annealing  produced  no  appreciable 
change  in  its  transmission.)  It  is  probable  that  the  increase  in  transmission  is  due  primarily  to 
the  removal  of  Te-rich  inclusions  or  microprecipitates  formed  during  growth.  This  explanation 
is  consistent  with  microscopic  observations  made  on  alloy  samples  which  contained  enough  MgTe 
to  be  transparent  to  visible  light.  The  as-grown  samples  contained  large  numbers  of  small  opaque 
areas  which  were  greatly  reduced  in  size  by  annealing. 

Transmission  measurements  between  2  and  15  pm  were  made  on  annealed  samples  at  room 
temperature  with  a  double-beam,  prism  spectrometer.  For  each  sample,  the  transmission  is 
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Fig.  111-5.  Absorption  coefficient  at  room  temperature  vs  photon  energy 
for  Cdj  ^Mg  Te  alloys. 


Fig.  111-6.  Energy  gap  at  room  temperature 
vs  composition  for  Cd^  ^Mg^Te  alloys. 
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almost  independent  of  wavelength  and  lies  between  50  and  65  percent.  These  results  indicate 

that  the  transmission  in  this  region  was  determined  primarily  by  reflection  losses. 

In  order  to  determine  the  energy  gaps  of  the  alloys,  absorption  coefficients  were  calculated 

from  the  transmission  data  in  the  vicinity  of  the  absorption  edge.  The  reflectivity  for  all  samples 

was  taken  to  be  0.24,  the  value  for  CdTe  at  its  absorption  edge.  The  results  are  shown  in  Fig.  Ill- 5 , 

where  the  absorption  coefficient  is  plotted  on  a  logarithmic  scale  as  a  function  of  photon  energy. 

The  results  for  CdTe  are  in  reasonably  good  agreement  with  those  reported  by  Marple.^  The 

absorption  edge  is  not  as  steep  as  Marple’s,  presumably  because  of  the  higher  purity  and  surface 

perfection  of  Marple's  sample  which  was  zone-refined  and  chemically  polished.  The  absorption 

edges  for  the  alloy  samples  are  considerably  less  steep  than  the  edge  for  CdTe.  This  difference 

could  easily  be  due  to  the  lower  purity  of  the  alloys.  Although  they  were  synthesized  from  the 

best  commercially  available  Mg,  this  material  is  much  less  pure  than  the  Cd  and  Te  used. 

An  arbitrary  but  consistent  method  was  used  to  evaluate  the  energy  gaps  of  the  alloys.  The 

energy  gap  for  each  sample  was  obtained  by  adding  0.  056  eV  to  the  photon  energy  for  which  the 

_  \ 

absorption  coefficient  of  the  sample  was  10  cm  .  The  correction  term  of  0.056eV  was  adopted 
in  order  to  yield  the  published  value,  1.50  eV,  for  the  energy  gap  of  CdTe  at  room  temperature. 

The  energy  gap  values  obtained  in  this  manner  are  plotted  against  mole-fraction  MgTe  in  Fig.  Ill  -  6 . 
The  data  are  well  represented  by  a  straight  line.  Extrapolation  to  x  =  1  gives  3.08  eV  as  the 
energy  gap  of  hypothetical  cubic  MgTe.  As  shown  in  Fig.  II I- 6 ,  energy  gap  values  calculated 

4 

from  the  absorption  edge  wavelengths  reported  by  Itoh  are  in  fairly  good  agreement  with  the 
present  values.  However,  the  present  values  are  not  consistent  with  Itoh’ s  conclusion  that  the 
absorption  edge  wavelength,  rather  than  the  energy  gap,  varies  linearly  with  MgTe  content.  Ex¬ 
trapolation  of  his  wavelength  vs  composition  line  gives  5.0eV  for  the  energy  gap  of  cubic  MgTe. 

7 

This  value  seems  much  too  high,  since  published  optical  absorption  data  for  thin  films  of  hex¬ 
agonal  MgTe  indicate  that  the  energy  gap  of  this  phase  is  less  than  3.6eV.  At  this  energy,  the 

5  -1 

absorption  coefficient  has  already  risen  to  about  10  cm  ,  and  there  is  a  peak  in  the  curve  of 


absorption  coefficient  vs  energy. 
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C.  PHOTOLUMINESCENCE  DUE  TO  OXYGEN  IN  ZnTe  ALLOYS 

g 

Thomas  and  co-workers  were  the  first  to  observe  mirror  luminescence  and  absorption 

9 

bands  in  ZnTe  due  to  a  center  lying  about  0.4  eV  below  the  conduction  band.  They  proposed 
that  this  center  is  an  isoelectronic  trap  resulting  from  the  substitution  of  an  oxygen  atom  on  a 

o 

tellurium  site.  At  low  temperatures,  the  luminescence  band  peaks  close  to  6500  A,  and  the  ab- 

o 

sorption  band  peaks  near  5900  A.  Both  bands  exhibit  marked  structure  due  to  phonon  cooperation. 

According  to  the  proposed  model,  the  oxygen  acts  as  an  acceptor-like  electron  trap,  since 
oxygen  is  more  electronegative  than  the  tellurium  it  replaces.  In  the  no-phonon  process,  the 
absorption  of  a  photon  of  energy  Eq  results  in  the  formation  of  a  bound  exciton,  which  is  com¬ 
posed  of  an  electron  strongly  bound  to  the  oxygen  center  and  a  hole  weakly  bound  to  the  electron. 
Luminescence  is  observed  when  the  exciton  is  annihilated.  In  the  no-phonon  case,  the  energy  of 
the  emitted  photon  is  again  Eq.  The  electron  trapping  energy  E^  is  equal  to  the  band  gap  energy 
Eg  minus  Eq.  For  ZnTe  at  4.2°K,  where  we  have  performed  our  luminescence  experiments, 

E  is  2.39eV,  E  is  1.986eV,  and  E,  is  therefore  0.40eV. 
go  t 
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Fig.  111-8.  Photoluminescence  energies  and  energy 
gap  for  ZnTe]_xSex  alloys  at  4. 2°K,  plotted 
against  mole-fraction  ZnSe. 


Fig.  111-7.  Spectra  of  photoluminescence  at  4.  2°K 
due  to  oxygen  in  ZnTej_xSex  alloys.  Wavelengths 
corresponding  to  no-phonon  energy  E  are  indi¬ 
cated  by  arrows. 


Fig.  IM-9.  Photoluminescence  energies  for 
ZnTe]  _XSX  alloys  at  4.  2°K,  plotted  against 
mole-fraction  ZnS. 
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At  the  beginning  of  our  investigation,  oxygen  had  not  been  conclusively  identified  as  the  cen¬ 
ter  responsible  for  the  mirror  luminescence  and  absorption,  although  Hopfield,  Thomas,  and 
9 

Lynch  had  reported  observing  the  most  intense  absorption  in  samples  grown  from  a  melt  doped 
with  ZnO.  Similarly,  we  have  found  that  the  absorption  coefficient  in  vapor-grown  crystals  in¬ 
creases  with  increasing  partial  pressure  of  oxygen  in  the  growth  ampoule.  Merz^°  has  recently 
presented  much  more  direct  evidence  for  the  oxygen  model  based  on  isotope  substitution  and 
Zeeman  experiments. 

In  this  investigation,  we  have  observed  the  variation  in  the  oxygen  luminescence  spectrum 
which  results  when  ZnTe  is  alloyed  with  ZnSe,  ZnS,  or  CdTe.  On  the  basis  of  the  oxygen  model, 
the  trapping  energy  is  expected  to  decrease  with  increasing  ZnSe  or  ZnS  content,  since  oxygen 
and  tellurium  differ  more  in  electronegativity  than  do  oxygen  and  selenium  or  oxygen  and  sulfur. 
On  the  other  hand,  is  expected  to  remain  essentially  constant  in  the  ZnTe-CdTe  alloys,  since 
changes  in  the  cation  lattice  should  not  strongly  affect  the  trapping  process.  These  predictions 
have  been  confirmed  by  the  data. 

Alloy  samples  were  prepared  by  annealing  weighed  mixtures  of  the  powdered  compounds  in 
sealed  quartz  ampoules.  In  order  to  enhance  the  luminescence,  excess  zinc  metal  was  added  to 
the  powder  mixtures,  and  the  ampoules  were  backfilled  with  oxygen  at  400  torrs.  The  samples 
were  annealed  at  1050°C  or  above  for  3  days  and  then  either  furnace-cooled  or  quenched  in  water. 
The  alloy  compositions  were  determined  from  x-ray  lattice  parameters  and  did  not  vary  from 
the  nominal  compositions  by  more  than  2  mole-percent. 

Figure  III— 7  shows  representative  photoluminescence  spectra  for  ZnTe-ZnSe  alloys  at  4.2°K 

o 

excited  by  3650  A  radiation.  The  luminescence  band  is  shifted  to  higher  energies  with  increasing 
ZnSe  content.  Addition  of  as  little  as  3  mole-percent  ZnSe  causes  the  phonon  structure  to  dis¬ 
appear.  This  makes  it  impossible  to  determine  Eq  for  the  alloys  from  the  position  of  a  distin¬ 
guishable  no-phonon  line.  On  the  basis  of  the  data  for  ZnTe,  E^  for  each  alloy  has  been  taken 
as  the  energy  corresponding  to  the  wavelength  on  the  high-energy  side  of  the  luminescence  band 
at  which  the  intensity  is  15  percent  of  the  maximum  intensity  for  the  band. 

In  Fig.  Ill- 8 ,  the  4.2°K  values  of  Eq  and  of  the  photon  energy  at  the  peak  of  the  luminescence 
band  are  plotted  against  x,  the  mole-fraction  of  ZnSe  in  the  ZnTe^^Se^  alloys.  Up  to  x  =  0.4, 
both  energies  increase  as  x  increases,  and  the  data  are  well  represented  by  two  parallel 
straight  lines.  This  indicates  that  the  shape  of  the  band  remains  essentially  the  same  over  this 
composition  range.  Near  x  =  0.4,  however,  there  appears  to  be  an  abrupt  change  in  the  character 
of  the  luminescence,  and  we  believe  that  the  luminescence  above  x  =  0.4  is  not  due  to  oxygen  but 
rather  to  some  unknown  impurity. 

Figure  III— 8  also  gives  the  values  of  the  energy  gap  E  for  the  ZnTe.  Se  alloys  at  4.2°K. 

g  l  x  x 

The  E  curve,  which  passes  through  a  minimum  at  x  =  0.4,  was  obtained  by  scaling  the  room- 

£  11 
temperature  values  determined  by  Larach,  Shrader,  and  Stocker.  It  is  seen  that  the  trapping 

energy  E^.,  which  is  the  difference  between  E^  and  Eq,  decreases  with  increasing  ZnSe  content, 

as  expected  for  the  isoelectronic  trap  model.  Extrapolation  of  the  Eq  line  indicates  that  the 

trapping  energy  becomes  zero  at  about  x  =  0.5.  Consequently,  there  is  no  luminescence  due  to 

oxygen  at  higher  concentrations  of  ZnSe.  For  these  concentrations,  the  luminescence  observed 

is  due  to  the  unknown  impurity  mentioned  above. 

Figure  III- 9  shows  the  oxygen  luminescence  data  for  ZnTe-ZnS  alloys  containing  up  to 
20  mole-percent  ZnS.  The  photon  energies  at  the  peak  of  the  luminescence  band  and  at  15-percent 
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intensity  increase  with  increasing  ZnS  content.  The  data  are  well  represented  by  two  parallel 
straight  lines,  whose  slope  is  somewhat  greater  than  the  slope  of  the  corresponding  lines  in  the 

ZnTe-ZnSe  system.  Although  energy  gap  data  for  the  ZnTe-ZnS  system  are  very  limited, 

1 1 

Larach,  Shrader,  and  Stocker  have  found 
that  1^^  initially  decreases  with  increasing  ZnS 
content.  Therefor^,  E^  decreases  with  in¬ 
creasing  ZnS  content,  as  expected  for  the  iso- 
electronic  trap  model. 

Figure  III-10  shows  the  luminescence  data 
for  the  ZnTe-CdTe  alloys.  In  addition  to  the 
oxygen  band,  a  second  band  at  somewhat  higher 
energy  was  observed  in  all  these  samples.  Be¬ 
cause  there  is  some  overlap  between  the  oxygen 
band  and  the  second  band  (which  is  presumably 
due  to  an  unknown  impurity  originating  in  the 
CdTe),  the  no-phonon  energy  for  the  oxygen 
band  cannot  be  determined  from  the  wavelength 
at  15-percent  intensity.  Therefore,  Fig.  Ill  - 1 0 
includes  only  the  photon  energy  at  the  peak  of 

the  oxygen  band,  the  energy  at  the  peak  of  the  second  band,  and  the  energy  gap  for  the  ZnTe- 

12  13 

CdTe  alloys,  which  is  reported  to  vary  linearly  or  almost  linearly  with  composition.  In 
contrast  to  the  results  for  the  alloys  with  ZnSe  or  ZnS,  the  energy  of  the  oxygen  luminescence 
band  decreases  with  increasing  CdTe  content.  The  data  are  well  represented  by  a  straight  line 
whose  slope  is  approximately  the  same  as  the  slope  of  the  energy  gap  Thus,  Et  remains  es¬ 
sentially  constant  over  the  composition  range  investigated,  as  expected  for  the  isoelectronic 

trap  model.  ^  T  i 

G.  W.  Iseler 

A. J.  Strauss 


Fig.  111-10.  Photoluminescence  energies  and  energy 
gap  for  Znj_xCdxTe  alloys  at  4. 2°K,  plotted 
against  mole-fraction  CdTe. 


D.  COMPOSITION-  AND  PRESSURE -INDUCED  STRUCTURAL  CHANGES 
IN  THE  Ba1_xSrxRu03  SYSTEM 

The  structural  chemistry  of  many  inorganic  compounds  with  the  general  formula  ABX^, 

where  A  and  B  are  cations  and  X  is  an  anion,  can  be  described  in  terms  of  hexagonal  and/or 

1 4 

cubic  close  packing  of  AX^  layers  with  the  B  cations  filling  the  X^  octahedra  thus  formed. 

When  all  layers  are  hexagonal  close  packed,  the  octahedra  share  faces  with  each  other  to  form 
chains  [Fig.  I II- 1 1  (a )] .  When  all  layers  are  cubic  close  packed,  the  octahedra  form  a  three- 
dimensional  array  by  sharing  only  corners,  giving  the  perovskite  structure  [Fig.  Ill—  11(b)]. 

If  structures  with  mixtures  of  hexagonal  and  cubic  close-packed  layers  are  allowed,  one 
might  expect  an  infinite  number  of  possible  structures.  However,  there  appear  to  be  only  three 
intermediate  stacking  sequences  in  actual  compounds  The  one  most  closely  resembling  cubic 
close  packing  has  a  6-layer  (6L)  sequence  cchcch  [Fig.  Ill- 1 1  (c )] ,  where  c  and  h  represent  cubic 
and  hexagonal  closed-packed  layers,  respectively  The  structure  having  equal  amounts  of  hex¬ 
agonal  and  cubic  close  packing  has  a  4-layer  (4L)  sequence  chch  [Fig.  Ill- 1 1  (d )] .  The  last  known 
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Fig.  II 1-1 1 .  Schematic  of  [1 10]  plane  of  hexagonal  cell  for  ABXg  structures  showing  face 
and  corner  sharing  by  octahedra  between  AX^  layers;  c  represents  c-parameter  of  crystal. 


intermediate  structure  has  predominantly  hexagonal  close  packing  with  the  9-layer  sequence 

15 

hhchhchhc  [Fig.  Ill-  11(e)].  Katz  and  Ward  have  discussed  in  detail  the  structural  relation  of 
all  these  types  plus  some  of  the  unknown  possibilities.  The  system  Ba1  xSrxRuO^  (Ref.  16)  is 
of  interest  since  one  is  able  to  study  the  effects  both  of  A  cation  size  and  of  pressure  on  the 
structural  transformations  from  the  9-layer  compound  (BaRuO^)  through  the  4-  and  6-layer 
phases  to  the  perovskite  (SrRuO^).  The  strontium-rich  members  are  also  interesting  since 
they  exhibit  spontaneous  magnetism. 

In  the  present  investigation,  ambient  pressure  material  was  prepared  by  the  solid  state 
reaction  of  BaCo^-SrCO^  mixtures  with  Ru  metal  at  temperatures  up  to  1100°C.  The  prepara¬ 
tion  of  high-pressure  phases  was  carried  out  in  cylindrical  gold  or  platinum  capsules  in  a  belt 
apparatus.  The  temperature  was  raised  after  application  of  pressure,  held  for  one-half  hour, 
and  then  rapidly  quenched  before  the  pressure  was  released.  Phase  and  structural  data  were 
obtained  by  analyzing  x-ray  diffraction  patterns  for  powdered  samples  measured  with  a  Norelco 
x-ray  diffractometer  using  monochromatic  CuKc*  radiation.  The  concentration  of  the  components 
in  two-phase  samples  was  estimated  from  the  relative  intensities  of  the  lines  in  the  diffraction 
patterns. 
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Fig.  111-12.  Pressure-vs-composition  phase  diagram  for  Ba^  ^Sr  RuO^  system  at  1000°C. 

The  results  of  this  investigation  are  summarized  in  the  pressure-vs-composition  plot  shown 
in  Fig.  Ill- 1 2  for  the  system  at  1000°C.  It  can  be  seen  that  the  addition  of  strontium  and  the  ap¬ 
plication  of  pressure  both  favor  structures  with  more  cubic  close  packing.  Both  effects  tend  to 
decrease  the  B-B  distances  and  therefore  to  make  the  face  sharing  of  octahedra  (associated  with 
hexagonal  close  packing)  less  stable.  In  the  6-layer  region  with  high  strontium  content,  distor¬ 
tion  of  the  normal  hexagonal  unit  cell  (a  s  5.6  A,  c  a  13.  9  A)  produces  diffraction  line  splitting. 
The  splitting  can  be  attributed  to  a  monoclinic  distortion  with  a  ^  b  and  y  ^  60°.  Even  in  the 
most  distorted  structures,  the  deviation  from  a  hexagonal  unit  cell  is  less  than  4  percent. 

J.  M.  Longo  D.  A.  Batson 
J.  A.  Kafalas  T.  W.  Hilton 

E.  DETERMINATION  OF  EXPERIMENTAL  FORM  FACTORS  ON  AN  ABSOLUTE 
SCALE  FROM  RELATIVE  X-RAY  INTENSITIES 

The  controversy  about  experimental  determination  of  atomic  form  factors  in  solids,  by 
means  of  x-ray  diffraction  intensity  measurements,  relates  only  to  results  on  an  absolute  scale. 
On  a  relative  scale,  all  the  available  experimental  results  agree.17  An  alternate  method  of 
evaluating  the  scaling  constant  would  thus  be  useful. 

The  form  factor  can  be  represented  by  the  parametric  expression 

f(s)  =  Yj  exP  f— 

j 
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where  s  =  |k|/47r  =  (sin0)/A,  with  0  the  Bragg  angle  and  \  the  wavelength,  and  the  a's  and  0's 

are  Gaussian  coefficients.  We  require  the  form  factor  to  satisfy  the  two  conditions  f(s  =  0)  =  Z, 

where  Z  is  the  atomic  number,  and  f(s  =  f°(s  >s c),  where  s c  is  the  reciprocal  lattice 

vector  value  at  which  the  outer  electrons’  contribution  to  the  scattering  is  negligible  and  f°(s) 

is  the  theoretical  free  atom  form  factor.  This  allows  a  least-squares  determination  of  the 

scaling  constant  as  well  as  the  Debye-Waller  temperature  parameter  B(T)  and  the  a’s  and  /3's. 

This  method  has  been  applied  to  our  relative  x-ray  intensity  data  on  Ni  metal  (Table  III  - 1 ) , 

1 8 

which  are  in  excellent  agreement  with  previous  results  collected  with  the  same  radiation 
(filtered  MoKcv  radiation).  Figure  III  - 1 3  shows  our  values  for  the  relative  structure  factor 
J(s)  =  (I(s)/mP(0)]  '  ,  where  I(s)  is  the  integrated  intensity  on  a  relative  scale,  m  is  the  plane 


TABLE  lll-l 

INTEGRATED  INTENSITIES 

OF  X-RAY  DIFFRACTION  LINES  FOR  Ni 

sin  9 

X 

hkl 

l(s)f 

l(s)  cal* 

§ 

P(s)9 

0. 2458 

111 

351 . 84 

349.  78 

351.01 

0.2838 

200 

169.  84 

169.06 

170.82 

0.4013 

220 

101. 13 

100.04 

99.60 

0.  4706 

311 

104. 65 

104. 43 

103.23 

0.4915 

222 

28.91 

28.  84 

28.  22 

0.5676 

400 

11.48 

11.42 

10.88 

0.6185 

331 

29.29 

29.74 

29.  26 

0.6346 

420 

26.03 

25.56 

25.70 

0.6951 

422 

16.56 

16.53 

16.  16 

0.  7373 

333,511 

16.21 

16.29 

15.88 

0.8027 

440 

3.80 

3.96 

3.  85 

0.8394 

531 

12.79 

12.80 

12.63 

0.8514 

442,600 

7.56 

7.38 

7.  58 

0.8974 

620 

4.45 

4.  57 

4.  44 

0.9305 

533 

3.95 

3.87 

4.  06 

0.9412 

622 

3.67 

3.68 

3.75 

0.9831 

444 

1.02 

1.02 

1.01 

1.0133 

551,711 

5.  59 

5.  44 

5.  58 

1.0232 

640 

2.73 

2.62 

2.75 

1.0618 

642 

4.61 

4. 62 

4.79 

1 . 0899 

553,731 

6.33 

6.41 

6.  60 

1. 1351 

800 

0.53 

0.48 

0.51 

1.1614 

733 

(1.89 

1.88 

1.92 

1.  1701 

644,820 

13.71 

3.69 

3.80 

1 . 2040 

660,822 

2.65 

2.69 

2.64 

1.2288 

555,751 

1 4.  17 

4.  34 

4.26 

1.2370 

662 

11.82 

1.80 

1.87 

1.2691 

840 

1.85 

1.83 

1.90 

1.2927 

753,911 

(5.79 

5.  73 

5.  95 

1 . 3005 

842 

13.91 

3.87 

4.02 

t  Observed  intensities.  The  lines  joined  by  a  brace  were  overlapping.  They 
have  been  counted  together  and  resolved  according  to  the  theoretical  ratio. 

t  Values  calculated  from  the  Gaussian  polynomi 

al  representing  the  data. 

§  Experimental  results  of  Ref.  18. 
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Fig.  111-13.  Relative  structure  factor  J(s)  as  function 
of  s  =  (sin  0)/X  for  Ni  metal. 


'  TABLE  II 1-2 

VALUES  OF  X-RAY  PARAMETERS 

FOR  Ni  BEFORE  AND  AFTER  REFINEMENT 

Initial 

Final 

c 

0.0020 

0.00166 

ai 

1 . 2640 

1 . 27705 

Pi 

0.0020 

0.11825 

a2 

13.2174 

13.2296 

P2 

3.  7504 

3.72131 

a3 

7.  1019 

7. 06452 

P3 

0.2671 

0. 24092 

°4 

4.1801 

4. 1 8757 

P4 

12.4514 

13.14709 

a5 

2.  2246 

2.24115 

P5 

60. 3462 

173.8109  7 

B(T) 

0.381 

0.3705 
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Fig.  111-14.  Atomic  form  factor  obtained  from  ex¬ 
perimental  data  for  Ni  metal  compared  with  theo¬ 
retical  free  atom  form  factor  calculated  by  Watson 
and  Freeman  (Ref.  20). 


Fig.  111-15.  Difference  (Ap)  between  charge  dis¬ 
tribution  obtained  from  experimental  data  for  Ni 
metal  and  theoretical  charge  distribution  for  free 
Ni  atom,  as  function  of  distance  (r). 


multiplicity,  and  P(9)  is  the  Lorentz  and  polarization  correction.  Figure  III- 13  also  shows  the 
Gaussian  line  best  representing  the  J  values.  Table  III  -  2  lists  the  parameters  before  and  after 

1 9 

least- squares  refinement.  The  initial  a  and  values  are  those  obtained  by  Cromer  and  Waber 
for  the  free  Ni  atom.  The  accuracy  of  the  fit  in  Fig.  Ill- 13  averages  2  percent,  the  same  as  our 
experimental  error  and  our  discrepancy  with  previous  results.  In  this  experiment,  the  major 
difficulty  is  in  making  the  background  determinations  with  sufficient  accuracy. 

8  2 

In  Fig.  Ill- 1 4,  the  theoretical  free  atom  form  factor  for  the  electronic  configuration  3d  4s 
(Ref.  20)  is  compared  with  the  best  Gaussian  representation  of  our  results  for  the  form  factor  of 
solid  Ni.  One  sees  that  there  is  a  systematic  deviation  which  reaches  a  maximum  of  3  percent 
for  s  ~  0.6.  This  is  small  but  probably  outside  the  experimental  error.  Fourier  inversion  of 
the  two  Gaussian  representations  yields  the  atomic  charge  distributions  for  the  free  atom  and 
for  the  solid.  The  difference  between  these  distributions  is  shown  in  Fig.  Ill- 1 5 ,  where  one  sees 

o 

that  the  atom  in  the  solid  has  more  electrons  than  the  free  atom  in  the  region  from  0.17  to  0.37  A, 

o 

and  less  between  0.37  A  and  infinity.  This  suggests  a  transfer  from  the  outer  levels  (4s)  to  the 

°  20 
first  inner  levels  (3d).  Since  0.32  A  is  the  position  of  the  3d  level  in  the  free  Ni  atom,  numerical 

integration  of  the  real  space  charge  density  suggests  a  transfer  of  0.3  electron,  a  value  smaller 
than  suggested  by  the  magnetic  moment.  However,  the  accuracy  of  this  determination  is  poor, 
and  the  result  should  be  considered  as  qualitative. 

It  appears  that  the  scaling  method  proposed  here  preserves  the  information  concerning  the 
outer  electron  distribution  which  is  lost  when  the  theoretical  free  atom  form  factor  is  used  in 
scaling  over  the  whole  experimental  range.  The  very  slight  difference  from  the  free  atom  dis¬ 
tribution  found  for  Ni  suggests  an  electronic  configuration  3d^*^4s^  ^ ,  where  6  =  0.3.  Further 
work  could  convert  this  qualitative  result  to  a  quantitative  value  having  a  known  accuracy. 


P.  M.  Raccah 
V.  E.  Henrich 
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IV.  PHYSICS  OF  SOLIDS 


A.  ELECTRONIC  BAND  STRUCTURE 

1.  Optical  Properties  of  Nickel-Copper  Alloys  from  0.2  to  11  eV 

Our  investigation  of  the  band  structure  of  nickel,  through  optical  studies,  has  led  us  to  the 

1  2 

nickel-copper  disordered  alloy  system.  Experimental  measurements  and  band  calculations 

indicate  an  energy  band  picture  for  copper  and  nickel,  around  the  L  point  in  the  Brillouin  zone, 

1 

as  shown  in  Fig.  IV- 1.  Also  given  here  are  the  low-energy  transitions  observed  in  nickel  and 
the  reversal,  between  copper  and  nickel,  of  the  L^  d-like  and  the  LI,  p-like  states.  It  seemed 
logical  to  test  the  band  model  for  nickel,  wrhich  is  less  well  established  than  that  of  copper,  by 
looking  at  the  optical  properties  of  the  alloy  sylstem  of  these  two  elements;  they  lie  next  to  each 
other  in  the  periodic  table,  copper  having  one  more  electron.  As  a  result,  we  find  that  the  alloy 
system  is  interesting  in  its  own  right. 

We  have  measured  the  reflectance  of  single  crystals  of  several  alloys  from  0.2  to  1 1  eV. 

Figure  IV-2  shows  the  reflectance  of  a  high-nickel  and  a  high-copper  alloy  along  with  pure  nickel 

and  pure  copper.  The  curves  for  the  pure  metals  are  virtually  the  same  as  those  reported  by 
1  3 

Ehrenreich,  et  al.,  and  Beaglehole.  The  0.4-eV  shoulder  in  nickel  is  absent  at  room  tempera¬ 
ture  in  the  61 -percent  nickel  alloy.  At  77  °K,  in  measurements  not  shown,  a  small  kink  in  the 
curve  for  the  61 -percent  nickel  alloy  does  seem  to  occur  at  0.3  eV;  this  may  be  related  to  ferro¬ 
magnetic  behavior.  We  have  previously  interpreted  this  structure  in  nickel  as  a  minority  spin 
band  transition  from  LI,  to  L^2  (see  Ref.  1). 


MINORITY  SPIN  BAND 


NICKEL 

Fig.  IV— 1 .  Energy  bands  of  nickel  and  copper  near 
L  point  in  Brillouin  zone. 


Fig.  IV-2.  Reflectance  at  room  temperature  of  nickel , 
copper,  and  two  alloys  of  these  metals. 
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In  the  high-copper  curves,  even  with  16-percent  nickel,  the  reflectance  starts  dropping  im¬ 
mediately  rather  than  staying  at  a  high  value.  In  the  three  samples  with  copper,  there  is  con¬ 
sistently  a  distinct  change  in  slope  near  2  eV,  characteristic  of  the  first  interband  transition  in 
copper  from  the  d-bands  to  the  Fermi  surface.  This  edge  does  not  shift  as  we  would  expect  if  an 
average  potential  or  rigid  band  model  applied  to  the  alloy  system.  In  either  of  these  models,  add¬ 
ing  nickel  to  copper  would  change  the  d-band  Fermi  surface  separation  and  shift  the  energy  of  the 

....  4 

first  interband  transition.  The  recent  limited  data  of  Schroder  and  Onengiit  on  alloys  with  as 
low  as  2-percent  nickel  also  show  a  drop  in  the  low-energy  reflectance  with  the  edge  remaining 
at  2  eV. 

At  higher  energies,  other  changes  occur.  The  61 -percent  nickel  alloy  shows  a  broad  peak 
occurring  between  2  and  5  eV,  with  a  distinct  narrowing  and  lowering  of  the  peak  at  5.2  eV.  A 
similar  but  not  so  distinctive  change  occurs  in  the  high-copper  samples. 

Our  results  are  related  to  band  models  by  using  a  Kramers-Kronig  analysis  to  compute  the 
imaginary  part  of  the  dielectric  constant  Figure  IV-3  shows  e^,  which  is  proportional  to  the 

joint  density  of  states  for  interband  transitions. 
The  suggestion  from  the  reflectance  data,  that 
the  2-eV  transition  in  copper  does  not  shift  as 
nickel  is  added,  is  confirmed  by  the  curves. 

The  nonshift  of  the  2-eV  copper  absorption  is 
clearly  evident  as  a  knee  in  the  16-percent  nickel 
alloy  spectrum  and  an  inflection  point  in  the  61- 
percent  nickel  alloy.  In  the  16-percent  nickel  al¬ 
loy,  a  rise  in  above  that  of  the  free  electron 
behavior  of  copper  is  obtained.  This  apparently 
comes  from  low-energy  transitions  to  the  Fermi 
surface  from  states  near  the  Fermi  surface. 

At  higher  energies,  the  detailed  picture  for 
the  alloys  is  less  clear  but  it  appears  that  other 
band  models,  such  as  the  "virtual  states"  model 


Fig.  IV-3.  Frequency  variation  of  imaginary  part 
of  dielectric  constant  e0 


5  6 

of  Friedel  and  Anderson,  and  the  "minimum 

polarity"  model  recently  discussed  by  Lang  and 
7 


Ehrenreich,  must  be  considered  in  order  to  ex- 
c2*  plain  the  data  below  4  eV.  The  "rigid  band"  model 

g 

is  inadequate.  Recent  optical  and  photoemission  results  of  Seib  and  Spicer  should  provide  com¬ 
plementary  information.  We  are  planning  further  measurements  on  high-copper  alloys  where  the 

change  in  optical  properties  is  most  drastic.  J.  Scouler 

J.  Feinleib 
J.  Hanus 


2.  Optical  Studies  in  Bismuth-Doped  Silicon 

9 

Measurements  (which  were  described  in  the  previous  Solid  State  Research  Report  )  have  been 
completed  of  the  optical  absorption  due  to  low-energy  transitions  of  isolated  bismuth  impurity 
atoms  in  silicon  at  low  temperatures  as  a  function  of  the  polarization  of  the  infrared  light  beams 
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g 

and  the  orientation  of  the  uniaxial  stress  applied  to  the  sample.  Uniaxial  stresses  up  to  6  X  10 
dynes/cm  were  applied  in  the  [001],  [110],  and  [111]  crystal  directions,  with  the  E  vector  of  the 
light  parallel  or  perpendicular  to  the  applied  stress.  An  example  of  the  data  is  given  in  Fig.  IV-4, 
showing  the  transmission,  in  arbitrary  units,  for  an  applied  stress  of  1.17  X  10  dynes/cm  in  a 
[001]  crystal  direction.  The  zero  stress  transition  energies  are  also  indicated  on  the  figure. 

All  measured  relative  intensities  and  transition  energies  as  a  function  of  strain  are  in  agree¬ 
ment  with  predictions  based  on  the  assumption  that  the  transitions  are  from  the  lsfA^)  ground  level 
(of  symmetry)  to  the  lsfT^)  level,  where  the  level  is  split  by  spin-orbit  coupling,  and  the 
form  of  the  strain  Hamiltonian  is  derived  from  effective  mass  and  deformation  potential  theory. 

W.  E.  Krag 
W.  H.  Kleiner 
H.  J.  Zeiger 


3.  Interband  Magnetoreflection  in  Gray  Tin 


Two  sets  of  interband  transitions  have  been  observed  in  the  magnetoreflection  of  gray  tin.  One 

set  has  been  seen  at  photon  energies  between  0.425  and  0.55  eV.  Points  from  these  transitions  ex- 

10 

trapolate  to  a  zero  magnetic  field  energy  gap  of  0.413  eV.  From  the  band  structure  of  gray  tin 
shown  in  Fig.  IV-5,  we  identify  these  transitions  as  taking  place  between  the  T-  valence  band  and 
the  Tg  conduction  band,  and  represent  them  by  arrow  A.  The  second  set  of  transitions  (arrow  B 
in  Fig.  IV-5)  has  been  observed  in  the  range  of  photon  energies  between  0.045  and  0.55  eV.  These 
take  place  from  the  T*  valence  band  to  the  conduction  band  across  the  zero  thermal  energy  gap. 


PHOTON  ENERGY  (eV) 


Fig.  IV-4.  Transmission  vs  photon  energy  of  bismuth- 
doped  silicon  for  uniaxial  appl ied stress  of  1.17  X10® 
dynes/cm^  in  [001]  crystal  direction  for  E  polarized 
parallel  and  perpendicular  to  applied  stress.  Zero  stress 
positions  of  absorptions  are  indicated  by  vertical  lines. 


Fig.  IV-5.  Band  structure  of  gray  tin  showing 
highest  valence  bands  and  lowest  conduction 
band  about  k  =  0.  Arrows  show  transitions  seen 
in  magnetoreflection. 
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14  -3 

For  this  experiment,  we  have  used  high-purity  gray  tin  which  has  nQ  ~  5X10  cm  ,  where 
nQ  is  the  density  of  donor  impurities  minus  the  density  of  acceptor  impurities.  This  impurity 
concentration  is  at  least  a  factor  of  10  lower  than  that  found  in  the  best  commercially  available 
tin.  The  improveme.nl  in  purity  is  accomplished  with  a  solvent  extraction  process  using  the  tin- 
mercury  alloy  system.  In  our  experiment,  light  is  reflected  off  a  natural  facet  of  a  crystal,  and 
most  of  our  measurements  have  been  made  on  a  crystal  with  a  (100)  face.  We  have  worked  in 
the  Faraday  configuration,  with  the  direction  of  light  propagation  parallel  to  the  magnetic  field  H; 
therefore,  our  data  are  characteristic  of  II  perpendicular  to  the  (100)  plane. 

The  magnetoreflection  results  for  gray  tin  are  very  similar  to  those  for  HgTe  (see  Ref.  11) 
but  are  more  interesting  for  a  comparison  with  theory  for  several  reasons.  First,  as  mentioned 
above,  they  are  characteristic  of  a  single,  high-symmetry  direction;  such  results  could  not  be 
achieved  in  HgTe  because  of  a  materials  limitation.  This  is  especially  important  for  fitting  the 
T g  Fg  transitions.  Second,  the  warping  of  the  F^  bands,  a  parameter  which  is  needed  in  the 
theory,  has  been  measured  for  gray  tin.1^  Third,  the  complications  which  may  be  introduced  by 
the  inversion  asymmetry  of  zincblende  HgTe  are  absent  for  gray  tin,  which  has  inversion 
symmetry. 

A  matrix  Hamiltonian,  which  includes  the  full  interaction  between  F  *  T  ",  and  T  *  bands 

2  ((o 

and  includes  higher  band  interactions  to  order  k  ,  has  been  used  to  calculate  the  magnetic  en~ 

13 

ergy  levels.  A  shortcoming  of  this  theory  for  both  gray  tin  and  HgTe  is  that  the  interaction  be¬ 
tween  Tg  and  the  conduction  band  (using  a  mixture  of  single-  and  double-group  notations)  is 
included  only  to  order  k^.  Thus,  while  the  nonparabolicity  of  the  Tg  conduction  band  caused  by 
the  interaction  between  T 7  and  Fg  bands  is  included,  the  nonparabolicity  of  the  band  caused  by 
the  r.,-  band  is  not. 

^  _  + 

Figures  IV- 6  and  IV-7  show  the  transition  energies  vs  magnetic  field  for  the  r_  -*■  Fft  and 
+  +  f  o 

T g  -*■  T g  transitions,  respectively.  The  solid  curves  are  theoretical  and  give  a  good  fit  to  both 
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0  20  40  60  80  100 

MAGNETIC  FIELD  (kOe) 

Fig.  IV-7.  Tq  -*■  Tg  transition  energies  showing  experimental  points  and  theoretical 
curves  with  same  convention  as  in  Fig.  IV-6.  These  transitions  have  been  seen  up  to 
0.55  eV,  but  are  plotted  here  to  only  0.2  eV. 

sets  of  transitions.  The  data  of  Fig.  IV-6  were  taken  with  the  sample  at  ~1.5°K,  and  those  of 

Fig.  IV-7  with  the  sample  at  ~20°K.  We  have  measured  the  transitions  of  Fig.  IV-6  at  20°  and 

~85°K  and  cannot  detect  any  shift  of  the  resonant  energies.  This  is  in  marked  contrast  to  HgTe 

14 

which  shows  large  temperature-induced  shifts. 

The  systematic  deviation  at  high  photon  energies  for  both  sets  of  transitions,  with  the  experi¬ 
mental  points  lying  below  the  theoretical  curves,  is  caused  by  the  nonparabolicity  of  the  con¬ 
duction  band  induced  by  the  T.  q  conduction  band,  as  mentioned  above.  This  effect  has  been  dis- 
11  1  ^ 

cussed  previously.  One  interesting  feature  of  the  low-energy  transitions  is  that  they  give  the 
valence  band  mass  directly.  This  is  evident  from  an  examination  of  the  Landau  levels  for 
these  transitions,  shown  in  Fig.  IV-8.  The  spin  splittings  of  T  *  valence  and  conduction  bands  are 
nearly  equal,  which  means  the  two  cr  (An  =  +1)  transitions  from  the  n  valence  sub-band  to  the 
n  +  1  conduction  sub-band  are  nearly  degenerate  in  energy.  The  same  is  true  of  the  a  (An  =  —  1) 
transitions.  Thus,  the  spacing  in  Fig.  IV-7  between  a  transitions,  or  between  a+  transitions,  is 
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Fig.  IV-8.  Lowest  Landau  levels  for  Tq  bands 
at  70  kG.  Arrows  show  optical  transitions  to 
n  =  2  levels  of  conduction  band.  Difference 
between  a  and  a  transition  energies  to  a 
given  conduction  band  level  is  direct  measure 
of  valence  band  level  spacing. 


the  usual  reduced  mass  cyclotron  energy.  However,  the  spacing  between  u  and  cT  transitions 
to  the  same  conduction  sub-band  is  twice  the  valence  band  cyclotron  energy. 

The  lowest  <7+  transition  in  Fig.  IV-7  is  at  too  high  an  energy  to  be  part  of  the  Landau  level 
transitions  shown  in  Fig.  IV-8.  Our  suspicion  is  that  it  is  a  transition  from  an  impurity  ground 
state,  associated  with  the  n  =  0  conduction  band  Landau  level  to  an  excited  impurity  state  asso- 

r 

ciated  with,  and  lying  below,  the  n  =  2  conduction  sub-band.  The  observation  of  impurity  transi- 

15 

tions  similar  to  this  in  p-type  InSb  and  germanium  makes  this  plausible  from  an  experimental 
point  of  view.  Evidence  has  been  presented  recently  for  discrete  impurity  levels  in  gray  tin  un¬ 
der  uniaxial  stress  which  causes  the  splitting  of  the  bands It  would  not  be  surprising  if  a 

similar  effect  were  caused  by  a  high  magnetic  field. 

11  + 

The  following  parameters  have  been  used  to  generate  the  theoretical  curves:  E(  )  — 

E(r^)  =  0.413  eV;  E(Tg)-E  (T+)  =  0.80  eV;  Ep  =  24.5  eV;  y^  =  4.42;  yz  =-1.35;  y3  =  1.35;  and 

k  -  —  3.0.  The  band  edge  masses  corresponding  to  this  set  of  parameters,  in  units  of  the  free 

electron  mass,  are  m  (T^)=  0.028,  m  ( T^)  =  — 0. 195,  m( )  =  —0. 057,  and  m( T*)  =  —  0.  051. 

In  addition  to  the  work  on  the  lowest  a+  transition  (Tg  -*•  T g  ),  we  are  hoping  to  obtain  more 

information  in  the  near  future  about  the  conduction  band,  from  its  nonparabolic  effect  on  the 

Tq  conduction  band,  and  to  make  a  detailed  comparison  of  our  magnetoreflectivity  results  with 

those  from  other  experiments,  especially  the  conduction  band  mass  determinations  from  the  os- 

cillatory  magnetoresistance.12- 17  s.  H.  Groves  A.  W.  Ewald 

C.  R.  Pidgeon  R.  J.  Wagner 
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4.  Inversion  Asymmetry  and  Warping-Induced  Interband 
Magneto -Optical  Transitions  in  InSb 

We  have  completed  a  study  of  "extra”  transitions  in  the  interband  magnetoreflection  and 
magneto-electroreflection  of  InSb,  caused  by  the  warping  and  by  the  linear-k  splitting  of  the  va¬ 
lence  band.  An  unambiguous  assignment  of  the  transitions  is  made  from  a  study  of  the  anisot¬ 
ropy  of  the  spectra,  using  left-  and  right-circularly  polarized  light,  with  the  magnetic  field 
along  [100],  [110],  [111],  and  [211]  directions  of  the  crystal.  From  the  relative  strength  and 
energies  of  these  transitions,  we  are  able  to  determine  the  size  of  the  warping  and  linear-k 
terms  in  the  effective  mass  Hamiltonian.  We  find  the  warping  parameter  (y^  — y =  1.0,  and 
the  linear-k  parameter  C  =  4.2eV-cm.  The  error  in  the  determination  of  (y^  —  y^)  is  about  15 
percent,  and  that  in  the  determination  of  C  is  about  30  percent.  The  result  for  (y~  —  y?)  is  in 

18  j  c. 

good  agreement  with  the  cyclotron  resonance  work  of  Bagguley,  et  al.,  but  not  with  the  more 

19 

recent  work  of  Tohver  and  Ascarelli.  The  value  of  C  is  about  three  times  smaller  than  sug¬ 
gested  by  Robinson,^  but  is  in  good  agreement  with  the  theoretical  estimate  of  Kane.  * 

22 

A  short  report  of  this  work  has  been  published  in  the  Physical  Review  Letters.  A  longer 
version,  with  full  account  of  the  experimental  and  theoretical  details,  is  being  prepared  currently. 

S.  H.  Groves 
C.  R.  Pidgeon 


5.  Determination  of  Debye  Temperatures  from  Thermal 
Modulation  of  X-Ray  Diffraction 

The  usual  method  of  measuring  the  Debye  temperature  9^,  in  terms  of  x-ray  diffraction  in¬ 
tensities,  consists  of  measuring  the  difference  between  the  Debye-Waller  factors  at  two  fairly 
23 

distant  temperatures.  This  method  precludes  the  possibility  of  making  detailed  studies  of  0^ 
in  any  small  temperature  range,  and  thus  is  not  useful  for  studying  critical  transitions.  We  have 
measured  the  temperature  derivative  of  the  Debye-Waller  factor  at  one  temperature  directly  by 
synchronous  detection  of  the  diffracted  x-ray  beam  from  a  thermally  modulated  sample.  As  the 
temperature  of  a  crystal  used  as  sample  in  an  x-ray  diffractometer  is  raised,  thermal  expansion 
of  the  lattice  shifts  the  Bragg  peaks  toward  lower  angle.  For  a  cubic  crystal,  this  can  be  writ¬ 
ten  as 

d(20o>  _  Xs /  h2  +  k2  +  i2  ,1  da, 

dT  ~  a  cos©  'a  dT '  1  ' 

o 


where  0q  is  the  Bragg  angle,  T  is  the  temperature,  X  is  the  x-ray  wavelength,  a  is  the  cell 
edge,  and  (hki)  is  the  diffracting  planes.  At  the  same  time,  the  increased  atomic  motion  changes 
the  intensities  of  the  peaks,  in  general  reducing  them.  For  a  simple  lattice,  this  is  taken  into 
account  by  the  Debye- Waller  factor  B(T),  which  multiplies  the  T  =  0°  intensity  in  the  form 

I(T)  =  1(0)  e"2B(T)  S  (2) 

where  S  =  sin0/x.  Thus, 


dl_  _  ofdB(T) 
dT  1  dT 


]  S2I(T) 


(3) 
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Fig.  IV-9.  Experimental  setup  for  thermal  modulation. 

If  the  sample  temperature  is  modulated  and  the  diffractometer  angle  is  held  constant,  the 
thermal  expansion  appears  in  the  AC  component  of  the  diffracted  beam  as  the  angular  derivative 
of  the  lineshape,  while  the  intensity  change  gives  a  signal  proportional  to  the  lineshape.  Anal¬ 
ysis  of  the  resulting  signal  yields  the  relative  amounts  of  these  two  components.  Then,  if  (l/a) 
(da/dT)  is  known,  dB/dT  can  be  determined. 

By  using  chopped  radiation  from  a  projection  lamp  focused  on  the  sample  (Fig.  IV-9)  to  mod¬ 
ulate  the  temperature  (modulation  amplitude  —  0.01°  to  0.1°C),  signals  of  the  form  shown  in 

Fig.  IV- 10  have  been  obtained  with  single-crystal 
NaCl.  Most  of  the  signal  is  due  to  thermal  expan¬ 
sion,  but  careful  examination  shows  the  presence  of 
a  small  contribution  from  the  changing  phonon  den¬ 
sity.  In  practice,  only  the  region  near  the  zero  of 
signal  is  measured,  since  the  location  of  the  zero 
uniquely  determines  the  amount  of  each  component 
present.  Several  runs  have  been  made  with  both 
CuKo  and  MoKo  radiation,  and  the  results  are  pre¬ 
sented  in  Table  IV- 1. 

When  the  theory  is  extended  to  a  diatomic  crys¬ 
tal,  an  average  Debye-Waller  factor  can  be  defined, 
and  this  is  the  quantity  whose  derivative  is  given  in 
Table  IV- 1.  Within  the  limits  of  the  Debye  theory, 
which  predicts  dB/dT  >  0  always,  the  usual  Debye 
temperature  can  be  determined  directly  from  dB/dT. 
This  has  been  done  for  the  CuKa  result,  and  there 
is  good  agreement  with  the  values  of  0D  =  290°K  determined  by  other  methods. 


| 3-SS-6861 | 


I _ I _ I _ I 

30.0  29  5  29  0  28.5 

DEGREES  26 

Fig.  IV— 1 0.  Thermal  modulation  signal  vs  20. 
Kaj  and  Ka2  were  partially  resolved  in  x-ray 
lineshape. 
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TABLE  IV-1 

dB/dT  COMPUTED  FROM  THERMAL  MODULATION  DATA 

(hki) 

Radiation 

dB/dT  (A/°C)  X  103 

od(°k) 

(400) 

(400) 

(600) 

(800) 

(10,0,0) 

CuKa 

MoKa 

MoKa 

MoKa 

MoKa 

+  5.11  ±1.2 

+  0.107  ±  1.24 

-0.76  ±  1.26 

— 1.79  +  0.70 

-2.34  +  0.76 

+  26 

281 

-31 

Under  any  reasonable  assumptions,  dB/dT  should  be  independent  of  (hki)  and  the  x-ray 

wavelength,  and  we  should  obtain  the  same  value  in  all  our  measurements.  We  believe  that  the 
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strange  results  obtained  with  MoKa  radiation  are  due  to  extinction.  CuKcv  radiation  penetrates 
only  about  60  pm  into  NaCl,  and  any  surface  damage  probably  extends  this  far  into  the  crystal, 
resulting  in  a  mosaic  structure  with  little  extinction.  However,  MoKa  radiation  extends  about 
560  pm  into  NaCl,  where  surface  damage  should  be  absent  and  extinction  could  be  expected  to 
occur.  Since  the  amount  of  extinction  depends  on  the  state  of  perfection  of  the  lattice,  we  believe 
that  the  thermal  gradients  set  up  in  the  sample  cause  periodic  variations  in  strain  which  modu¬ 
late  the  strength  of  the  extinction  and  hence  the  intensity  of  the  diffracted  beam.  This  effect  has 
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been  observed  directly  in  quartz,  and  its  size  could  easily  be  orders  of  magnitude  larger  than 
the  Debye-Waller  factor  contribution.  Tests  are  currently  under  way  to  determine  the  amount  of 
extinction  quantitatively. 

Several  criteria  affect  the  accuracy  obtainable  from  this  type  of  measurement.  Clearly,  ex¬ 
tinction  effects  must  be  eliminated.  Sandblasting  of  the  crystal  surface  and  use  of  very  soft 
x-rays  which  do  not  penetrate  deeply  into  the  sample  should  overcome  this  difficulty  in  most  ma¬ 
terials.  Also,  materials  with  lower  thermal  conductivities  will  have  larger  temperature  oscilla¬ 
tions  and  hence  larger  signals.  The  use  of  faster  detector  and  electronics  would  permit  use  of 
a  larger  photon  flux  (we  are  limited  to  30,000  photons  per  second).  Digital  rather  than  analog 
measurement  of  the  signal  would  both  simplify  data  analysis  and  increase  accuracy. 

With  the  improvements  mentioned  above,  this  technique  could  be  improved  to  the  point  of 
surpassing  existing  DC  methods.  In  addition,  the  ability  to  measure  0^  at  only  one  temperature 
gives  this  method  excellent  resolution  in  temperature,  making  it  valuable  in  the  study  of  phonon 

densities  around  critical  transitions.  Tr  ^  TT  .  , 

V.  E.  Henrich 

P.  M.  Raccah 


6.  Microwave  Acoustic  Amplification  in  n-InSb  at  9  GHz 

Acoustic  amplification  was  observed  on  applying  a  pulsed  electric  field,  up  to  40V/cm,  to  a 
rod-shaped  sample  of  n-InSb  at  4.2°K.  The  electric  field,  which  was  on  during  the  transit  time 
of  the  acoustic  pulse  down  the  sample  (about  2^sec),  produced  a  maximum  observed  change  in 
attenuation  of  —  41  db/cm.  Typical  results  are  shown  in  Fig.  IV-11  for  an  applied  field  of  40V/cm. 
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(  1  )  E  =  0 


(  2  )  E  =  40  V/cm 


A  — CAVITY  EXCITATION  PULSE 
B-  VIDEO  DETECTOR  OUTPUT 


Fig.  IV- 1 1 .  Output  of  receiver  video  detector  with  and  without  electric  field  applied. 

The  observed  amplification  was  enhanced  by  a  small  magnetic  field  (typically  less  than  50  gauss) 

applied  perpendicular  to  the  current  flow  which  presumably  eliminated  the  effect  of  injection  by 

the  end  contacts.  This  small  field  is  sufficient  to  move  any  injected  plasma  away  from  the  center 
2  6 

of  the  rod  (Suhl  effect)  and  out  of  the  1-mm-diameter  path  of  the  acoustic  pulse.  Since  o>cr  «  1 
for  B  =  50  gauss,  this  small  magnetic  field  was  neglected  in  our  theoretical  calculations.  One 
end  of  the  rod-shaped  sample  was  coated  with  aluminum  on  which  a  CdS  transducer  was  evapo¬ 
rated.  This  end  extended  into  a  9.05-GHz  coaxial  re-entrant  cavity  which  was  excited  by  a  5-W 

microwave  pulse,  y  to  1  psec  wide.  The  receiver  was  the  highly  sensitive  interferometer  de- 

27  13—3 

scribed  previously.  With  a  sample  electron  concentration  of  n  =  5  X  10  cm  ,  the  mobility 

2  8 

at  4.2 °K  is  strongly  electric  field  dependent  below  E  w  0.1  V/cm,  approaching  the  nitrogen  mo¬ 
bility  at  high  electric  fields  E  >  lOV/em.  In  addition,  the  electron  temperature  is  affected  by 
the  electric  field  for  fields  of  the  order  of  lOV/cm.  These  two  effects  complicate  the  theoretical 
interpretation  of  the  data.  Since  most  of  the  data  were  taken  with  electric  fields  above  E  =  1  V/cm, 
we  have  used  a  value  of  mobility  (i.e.,  qi,  where  q  is  the  sound  wave  vector,  and  l  is  the  elec¬ 
tron  mean  free  path)  in  the  theoretical  calculations  appropriate  to  the  higher  electric  fields. 

It  should  be  noted  that  the  zero-field  attenuation  decreases  about  5db/cm  on  increasing  the 
temperature  to  T  =  16°K  (see  Ref.  29),  where  the  attenuation  reaches  a  minimum.  Hence,  if  the 
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pulsed  electric  field  were  simply  heating  the  sample,  an  apparent  amplification  would  result. 

This  effect  was  minimized  by  bonding  one  contact  end  of  the  sample  to  the  copper  cavity  with 
indium  and  by  using  low  duty  cycle  (100-Hz)  short  electric  field  pulses.  Also,  in  order  to  moni¬ 
tor  the  temperature,  a  carbon  resistor  thermistor  was  placed  in  good  thermal  contact  with  the 
sample. 

Calculation  of  the  amplification  due  to  the  free  electrons  requires  essentially  only  a  knowl¬ 
edge  of  the  longitudinal  dielectric  constant  € .  .  The  lattice  dispersion  relation  for  a  piezoelectric 

30  * 

semiconductor  can  be  written 


CO 


=  qv 


so 


(4) 


where  co  and  v  are  the  frequency  and  sound  velocity  in  the  insulating  crystal  if  the  piezoelectric 
so  2  2  2 

coupling  is  vanishingly  small  (K^  «  1).  and  are  the  dimensionless  piezoelectric  and  de¬ 

formation  potential  coupling  parameters,  k€q  is  the  lattice  contribution  to  the  dielectric  constant, 
is  the  total  longitudinal  dielectric  constant  e ^  =  k€  —  cr/ioj,  and  cr  is  the  electrical  conduc¬ 
tivity.  Since  and  are  small,  the  attenuation  constant  ot  can  be  written 


“  =  2qi  =  ~<Kp  +  Kd2)  qIm(~e^)  (5) 

where  q.  is  the  imaginary  part  of  the  wave  vector,  and  Im  refers  to  the  imaginary  part.  On  sub- 

1  O  O  O  A  q  a 

stituting  the  value  ofK  =Kp+K^=3.6X10  obtained  previously, 


a 


2.26  X  io2Im(— db/cm 


(6) 


The  longitudinal  dielectric  constant  was  calculated  assuming  a  shifted  Maxwellian  electron 

distribution  with  the  effect  of  the  DC  electric  field  on  the  AC  response  included  through  the  use 

3 1 

of  a  complex  temperature.  We  note  that  the  classical  Boltzmann  equation  results  for  the  am¬ 
plification  in  our  experiments  differ  from  the  quantum  mechanical  density  matrix  results  by  less 
than  0.0025  percent,  despite  the  size  of  the  quantity  ftq/2mVo  ~  1.  Also,  note  that,  while  the 
complex  temperature  results  from  an  exact  solution  of  the  AC  Boltzmann  equation,  the  density 
matrix  equation  does  not  yield  this  simple  result.  As  in  the  classical  calculation,  the  quantum 
mechanical  solution  for  cot  »  1  includes  the  electric  field  in  the  AC  conductivity  exactly  through 
the  use  of  a  drifted  distribution  and  a  complex  temperature.  However,  when  diffusion  current 

is  present  (i.e.,  cot  <  1),  the  quantum  mechanical  solution  does  not  include  the  DC  electric  field 
32 

in  this  simple  manner.  Experimental  data  on  the  change  in  attenuation  with  electric  field  are 

shown  in  Fig.  IV- 12.  We  are  unable  to  measure  the  absolute  attenuation  at  zero  electric  field; 

however,  the  calculated  value  lies  between  6  and  10  db/cm  for  our  measured  low-field  mobilities 

(p  <  100,000  cmVv  sec).  The  initial  sharp  decrease  in  the  attenuation  is  probably  the  result  of 

the  rapid  increase  of  mobility  with  electric  field  in  this  range. 

The  results  of  the  theoretical  calculations  are  shown  in  Fig.IV-13.  The  value  of  ql  =  2 

corresponds  to  a  mobility  p  »  200,000  cm  /V  sec  which  is  nearly  the  experimentally  observed 

mobility  at  this  temperature  for  fields  >1  V/cm.  The  mobility  and  concentration  at  T  =  77 °K 
2  1 3  —  3 

were  p  =  630,000  cm  /V  sec  and  n  =  5  X  10  cm  ,  respectively. 
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Fig .  IV— 12 .  Experimental  change  in  attenuation  with  Fig.  IV-13.  Theoretical  dependence  of 

electric  field.  amplification  on  electric  field. 


The  considerable  discrepancy  between  the  theoretical  curve  and  experimental  data  may  re¬ 
sult  from  a  variety  of  effects,  among  which  are  injection  at  the  contacts  and  the  electric  field 
dependence  of  the  electron  mobility  and  temperature.  The  latter  has  not  been  considered  and 
this  may  introduce  a  considerable  error.  K  w  at.ii 


B.  MAGNETISM 

1.  Magnetic  Properties  of  MnAs  in  the  High-Pressure  Phase 

We  have  continued  our  study  of  the  magnetic  properties  of  the  high-pressure  (orthorhombic) 

phase  of  MnAs.  Our  measurements  of  the  magnetic  moment  vs  temperature  in  a  magnetic  field 
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of  lOkOe,  as  reported  in  the  previous  Solid  State  Research  Report,  have  been  extended  to  a 
pressure  of  llkbars.  In  addition,  we  have  now  studied  the  magnetic  moment  vs  temperature  in 
various  applied  fields  throughout  the  pressure  range.  A  set  of  such  measurements  (taken  at  a 
pressure  of  9  kbars),  as  given  in  Fig.  IV- 14,  illustrates  the  flattening  of  the  magnetization  curves 
in  an  intermediate  temperature  region,  as  well  as  the  onset  of  a  low-temperature  transition.  Re- 
expressing  the  data  in  terms  of  isotherms,  as  shown  in  Fig.  IV-15,  enables  us  to  define  the  mag¬ 
netic  state  of  the  sample  in  pressure-temperature  (P-T)  space.  Thus,  Fig.  IV-15  illustrates  the 
apparent  ferromagnetism  in  the  intermediate  temperature  region  as  well  as  the  metamagnetic  be¬ 
havior  of  orthorhombic  MnAs  at  both  higher  [Meta  (1)1  and  lower  [Meta  (2)]  temperatures.  These 
two  regions  are  distinguished  from  each  other  by  the  fact  that  in  state  Meta  (1),  dH^/dT  <  0, 
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Fig.  IV— 14.  Variation  of  magnetic  moment 
of  orthorhombic  MnAs  at  9  kbars  with  tem¬ 
perature,  for  various  applied  magnetic  fields. 


Fig.  IV- 1 5 .  Magnetic  moment  of  orthorhombic 
MnAs  at  9  kbars  as  function  of  applied  field  at 
various  temperatures,  indicating  both  ferromag¬ 
netic  and  metamagnetic  behavior. 
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Fig.  IV-16.  Magnetic  phase  diagram  for  MnAs. 


Fig.  IV-17.  Inverse  susceptibility  vs  temperature 
curves  of  orthorhombic  MnAs  at  indicated  pressures. 


TABLE  IV-2 

VALUES  OF  EFFECTIVE  MOMENT  AND  CURIE  TEMPERATURE 

OF  ORTHORHOMBIC  MnAs  AT  VARIOUS  PRESSURES, 

AS  DETERMINED  FROM  PARAMAGNETIC  SUSCEPTIBILITY 

Pressure 

Meff 

0p 

(kbars) 

(Bohr  magnetons) 

(°K) 

5 

1.91 

220.5 

6 

1.87 

228 

7 

1.83 

238 

8 

2.02 

242.5 

9 

2.02 

248.5 

10 

2.06 

256 
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while  in  Meta  (2),  dHcr/dT  >  0;  here,  ^cr  is  the  critical  transition  field  for  departure  from  anti¬ 
ferromagnetic  behavior.  Curves  similar  to  those  shown  in  Figs.  IV- 15  and  IV- 1 6  were  obtained 
for  various  pressures  and  temperatures.  These  enabled  us  to  define  the  region  in  P-T  space  in 
which  the  high-pressure  phase  is  ferromagnetic  or  one  of  the  metamagnetic  states.  These  regions 
are  delineated  in  Fig.  IV- 16. 

34 

In  addition.  Fig.  IV- 16  shows  the  hexagonal-rhombohedral  phase  transitions,  the  Neel  tem- 
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peratures  (T^t),  and  the  paramagnetic  Curie  points  9p.  These  latter  points  are  based  on  para¬ 
magnetic  susceptibility  data  taken  just  above  T^  in  the  pressure  range  from  4  to  lOkbars.  At 
these  temperatures  and  pressures,  MnAs  is  still  in  the  B31  phase.  The  resultant  curves  of  in¬ 
verse  susceptibility  in  this  region  follow  a  Curie-Weiss  relationship,  as  seen  in  Fig.IV-17,  and 
the  corresponding  values  of  \± ^  and  0p  are  given  in  Table  IV-2.  The  variation  of  observed 
over  this  pressure  range  is  small  and  falls  within  the  experimental  uncertainties.  However,  the 
value  of  approximately  2  pp  is  far  below  the  value  =  4.56  pp  which  was  observed  in  the  hexag¬ 
onal  phase  of  MnAs  35  This  indicates  that  virtually  the  entire  change  in  occurs  discontinu- 

ously  at  the  first-order  hexagonal-orthorhombic  transition,  and  directly  confirms  the  hypothesis 

3  6 

of  Goodenough,  Ridgley,  and  Newman  that  this  transition  is  accompanied  by  a  change  from  a 
high-  to  a  low-spin  manganese  configuration. 

The  paramagnetic  moment  value  is  consistent  with  the  saturation  moment  value  of  slightly 
over  1  |ig  obtained  in  the  ferromagnetic  state  at  llkbars.  In  addition,  the  positive  d0p/dP  indi¬ 
cates  a  net  increase  in  ferromagnetic  interactions  with  increasing  pressure,  which  is  also  con¬ 
sistent  with  the  magnetic  behavior  observed  below  the  ordering  temperature. 

N.  Menyuk 

J.  A.  Kafalas 

K.  Dwight 

2.  Magnetic  Resonance  in  Spiral  Spin  Structure  Compounds 

In  an  effort  to  gain  some  insight  into  the  nature  of  the  normal  modes  of  a  conical  spiral  sys¬ 
tem,  we  have  studied  the  equilibrium  configurations  and  spin  wave  spectra  for  the  linear  chain 
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model  mentioned  in  the  previous  Solid  State  Research  Report.  In  this  model,  two  species  of 
spins,  S^  and  Sp,  are  coupled  by  nearest  and  next-nearest  neighbor  A-B  and  B-B  interactions. 

It  is  found  that  there  exist  ranges  of  exchange  interactions  for  which  a  conical  spiral  ground 
state  is  stable.  Spin  wave  spectra  have  also  been  computed  for  a  number  of  conical  spiral  cases. 
Figure  IV-18  shows  the  spin  wave  spectrum  for  the  particular  choice  (in  normalized  units) 

=  3.5,  J^p  =  1.0,  Jpp  =  2.0,  Jpp  =  1.0,  and  S^  =  1,  Sp  =  0.8.  These  parameters  lead  to 
a  ferrimagnetic  conical  spiral  structure,  with  cone  angle  ip ^  =  5.81°,  =  121.9°,  and  turn  angle 

74.8°. 

The  modes  of  particular  interest  in  a  microwave  resonance  experiment  are  those  at  k  =  kQ, 
and  k  =  0.  For  the  k  =  kQ  modes,  those  marked  1,  2,  and  3  in  Fig.  IV-18  are  exchange  modes, 
while  the  mode  marked  4  is  the  uniform  mode.  For  modes  1  and  4,  each  conical  spiral  precesses 
rigidly  about  the  cone  axis,  while  for  modes  2  and  3,  the  cones  distort  as  they  precess  about  the 
cone  axis  (see  Fig.  IV- 19).  Modes  2,  3,  and  4  all  have  net  oscillating  moment  in  the  plane  perpen¬ 
dicular  to  the  cone  axes,  while  mode  1  has  no  net  moment. 
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Vk 


Fig.  IV-18.  Spin  wave  spectrum  for  linear 
chain-conical  spiral  configuration. 


Fig.  IV-19.  Schematic  diagram  of  spin  wave  modes 
of  linear  chain-conical  spiral  configuration.  Arrows 
show  instantaneous  displacement  of  spins  in  plane 
perpendicular  to  spiral  axis. 


MODE  4  MODES  2,3 

UNIFORM  PRECESSION  EXCHANGE 


MODE  1 
EXCHANGE 


MODE  b  MODE  6 

BREATHING  TWISTING 

AND 

TWISTING 
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For  the  k  =  0  case,  mode  6  is  a  pure  "twisting"  mode,  with  the  two  cones  moving  rigidly,  and 
therefore  to  =  0.  Mode  5  is  a  "breathing-twisting"  mode  (see  Fig.  IV-19),  for  which  the  net  os¬ 
cillating  moment  in  the  plane  vanishes.  Each  cone  separately  has  an  oscillating  moment  along 
its  cone  axis,  but  the  two  oscillating  moments  exactly  cancel. 


J.  J.  Stickler 
H.  J.  Zeiger 


3.  Spherical  Model  as  the  Limit  of  Infinite  Spin  Dimensionality 

Consider  the  Hamiltonian 


/(r) 


=  -4  J 


N 

Z 

i,  j= 


v.  .S.(l/)  -  S(f) 
lJ  1  J 


(7) 


where  =  [o',  (j),  a7(j),  .  .  .  ,  a  (j)]  are  r-dimensional  vectors  of  magnitude  \Ti7.  Although  it  is 

perhaps  not  generally  appreciated,  it  seems  clear  that  Jv  reduces  to  the  S  =  l/2  Ising,  classi¬ 
cal  planar,  and  classical  Heisenberg  models  for  v  =  1,  2,  and  3,  respectively. 

Here  we  argue  that  in  the  limit  v  -*  the  free  energy  of  approaches  that  of  the  spheri¬ 
cal  model.  This  result  is  of  particular  current  interest,  besides  the  geometrical  interpretation 
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it  attaches  to  the  spherical  model,  because  of  recent  evidence  that  various  "critical  properties" 
of  are  monotonic  functions  of  v.  Hence,  the  critical  properties  of  the  fairly  realistic  but 
hopelessly  insoluble  Heisenberg  model  (three-dimensional  spins)  would  appear  to  be  bounded  on 
one  side  by  those  of  the  Ising  model  (one-dimensional  spins)  and  on  the  other  by  those  of  the 
spherical  model  (infinite-dimensional  spins).  Moreover,  the  spherical  model  has  in  the  past  been 
interpreted  as  a  soluble  approximation  to  the  Ising  model,  whereas,  in  fact,  the  spherical  model 
would  appear  to  be  a  much  better  approximation  to  the  more  "realistic"  Heisenberg  model. 

The  normalized  partition  function  corresponding  to  is 


Q^)(K)  =  z^'uo/z^o) 


(8) 


where 


r-i 

*0O 

dci1(l),  da2(l 

oo  d, 

-oo 

N 

V 

x  n 

<5 

v  ~  Z  ffn(i) 

j=i 

n  =  l 

(N) 


Xexp[K2  v  2^=1an(i)  crn(j)J 


1J 

and  K  =  j/2kT.  The  N  constraints  are  then  represented  as  N  delta  functions,  and  we  obtain 

-*0'  +i°° 


(9) 


CK  +i°o 
-i« 


•  f  dt1 . dtN  exp [yF( K,  (t  })] 

'J  rv  -  i  oo  J 


(10) 


where  F(K,  |f})  =  KS.t.  —  [l/2  lnA({L})],  and  where  A({f})  is  the  determinant  of  the  quadratic 

form  2 .  (t  .6 . .  —  v. .)  a  (i)  a  (i). 
ij  J  ij  ij  n  n  J 
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The  requirement  that  F(K,  {V})  be  stationary  leads  to  the  set  of  conditions 

K  =  T  (9^  t  ;  i=1.2 . N  (11) 

Next,  expand  F(K,  { t .} )  about  the  stationary  points  {t  }, 

3  s 

N  /  2  \ 

F(K,{t})=F(K,{ts})  +  ^r  l  aflr)  '  (12) 

J  ,  m  n 

m,  n  =  l 

We  then  factor  yF(K,  {tg})  out  of  the  integrand  of  Eq.  (10),  calling  the  remaining  integral  R.  Thus, 
from  Eq.  (10), 

lnZ^fK)  =  j  pN  ln[7r/K]  +  i/F(K,  {t})  +  In  R  (13) 

The  limiting  free  energy  is  related  to  Z^'\k)  by 

-M(0O)  h  lim  lntZjJj^dO/Z^fO)]  (14) 

v ,  N~*°° 

where  (3  =  (kT)  *.  From  Eq.(13),  we  obtain 

-^(oo)  ln2K  +  |j  F(K,  (t  })  +  lim  In  R  (15) 

v,  N—°° 

The  last  term  on  the  right-hand  side  of  Eq.  (15)  can  be  shown  to  be  zero.  Now,  F(K,  {tg})  = 

NKt  -  [l/2  lnA({t  })]  and  lnA({t  })  =  ln2N  +  In  |(l/2  t  6..) -(l/2  v..)|.  Thus,  Eq.  ( 15)  is  pre- 
s  s  s  S  13  1J 

cisely  the  spherical  model  free  energy. 

H.  E.  Stanley 

4.  Eigenvalue  Degeneracy  as  a  "Mathematical  Mechanism"  for  Phase  Transitions 

39 

Kac  has  noted  a  striking  analogy  between  the  Ising  model 

HT=-J  Yj  S.(1)  S.(1)  (16) 

<i.  i> 

(with  strong,  short-range  forces)  and  the  "Kac  model" 

Hk  =  -yJ  £  e"YlHlSi(1)S.(1)  (17) 

h  j 

(with  weak,  long-range  forces  for  y  small).  In  the  Ising  model,  the  summation  is  over  nearest- 
neighbor  pairs  of  sites  <i,  j>  ,  whereas,  in  the  Kac  model,  the  summation  is  over  all  sites  i  and 
j  in  the  lattice.  Kac  observed  that  (for  both  models)  the  onset  of  a  spontaneous  magnetization 
below  a  certain  transition  temperature  Tc  is  characterized,  mathematically,  by  the  asymptotic 
degeneracy  of  the  maximum  eigenvalue  of  a  certain  linear  operator. 

Since  the  "mathematical  mechanisms"  responsible  for  the  phase  transitions  in  a  strong,  short- 
range  force  model  and  in  a  weak,  long-range  force  model  are  identical  for  "one-dimensional" 
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(i.e.,  Ising)  spins,  it  seems  natural  to  inquire  whether  this  same  mathematical  mechanism  is 
responsible  for  the  phase  transitions  in  other  spin  dimensionalities.  To  this  end,  we  consider 
the  Hamiltonians  obtained  by  replacing  the  one-dimensional  scalar  spins  S'  1  in  Eqs.(l6)  and 
(17)  by  isotropically  interacting  dimensional  vector  spins  and  obtain 


K^U-J  £  s!v)  ■ 

<i.j>  1 


(18) 


s  .  1,9) 

b  j 

Thus  far,  we  have  observed  that,  for  Eq.  (18)  with  v  =  1,  2,  3,  and  °°,  the  onset  of  infinite- 
range  order  may  indeed  be  described  mathematically  by  the  collapse  of  the  entire  (discrete) 
eigenvalue  spectrum  of  an  appropriate  linear  operator  into  a  single  value.  In  particular,  the 
spin-spin  correlation  function  may  be  expressed  as  ,  where  is  the 

largest,  and  X^  the  next-largest,  eigenvalue.  Hence,  the  spontaneous  magnetization  M  = 
lim  •  S^)>  is  zero  so  long  as  <  X^\  but  it  becomes  nonzero  as  the  eigenvalues 

become  degenerate.  H  E.  Stanley 

5.  Finite  Temperature  Single -Slater  Determinant  Theory 

We  have  applied40  standard  temperature-dependent  Hartree-Fock  (HF)  theory  (without  any 
symmetry  restrictions)  to  a  system  of  N  infinitely  separated  hydrogen  atoms  in  zero  external 
fields.  We  find  that,  as  temperature  T  -*•  0,  it  leads  to  the  entropy  S  -*■  0  rather  than  to  the  cor¬ 
rect  value  S  -*■  Nk  In  2,  although  the  zero-temperature  energy  and  the  absence  of  spin  ordering 
are  correctly  predicted.  Another,  closely  related  failure  of  that  theory  is  that  as  T  -*  0,  the 
spin  susceptibility  x  ®  rather  than  <*>  as  in  the  (correct)  Curie  law.  This  serious  failure  was 
traced  to  the  use  of  a  one-electron  operator  for  the  effective  or  approximate  Hamiltonian  H 
[defined  in  Eq.  (24)  below].  Furthermore,  a  new  finite-temperature  HF  theory,  still  involving 
single  determinants  as  approximate  energy  eigenstates,  was  shown  to  give  the  correct  behavior 
for  the  system  of  isolated  hydrogen  atoms. 

We  describe  here  some  of  the  essential  details  of  this  new  generalization  of  HF  theory. 

41 

Basic  is  the  well-known  variational  principle 

F(p)  =  tr(H  -  pN)  p  +  kT  trp  lnp  >  FT  (20) 


where  the  true  free  energy,  appropriate  to  a  grand  canonical  ensemble,  is 
Ft  =  -kT  In  tr  exp[-/?(H  -  |iN)] 


(21) 


p  is  any  density  operator,  (3  =  l/kT,  and  H  is  the  Hamiltonian 


H 


v+  +  1  _  +  + 

St..c.  c.  +  tt-  2  v.  .  i  .c.  c.  c.c. 
ij  i  3  2  13, ki  1  3  f  k 


(22) 


c.+  are  the  usual  fermion  creation  operators  for  some  set  of  orthonormal,  one-electron  states 
(p  The  equality  of  Eq.(20)  is  attained  when  p  is  the  true  density  operator  p^,  where 
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PT  =  Z'1  exp[-/3(H  -  pN)]  .  (23) 

It  is  convenient  to  define  an  approximate  Hamiltonian  H  such  that 

p  =  exp[-/3(H  -  pN)]/tr  exp[-/3(H -pN)]  .  (24) 

Standard  T-dependent  HF  theory  is  derived  from  Eq.  (20)  as  follows.  First,  H  is  restricted 
to  be  a  one-electron  operator 

H  e  2  CjN.  (standard  HF  theory)  (25) 

where  N.  are  the  occupation  number  operators  for  a  set  of  one-electron  states  (called  "orbitals,” 

involving  both  space  and  spin,  in  general),  and  are  c-numbers.  The  0.  and  are  then  chosen 

to  minimize  F(p),  yielding  the  standard  T-dependent  HF  equations 

(! I>v  ti^j)  +  2k  [(i^k,  v^^k)  -  V^k^j)l  Nk  =  <4<5i;j  (26) 

where 

Nk  =  tr  p  Nk  =  l/exp[(3(ei  -  p)  +  1]  .  (27) 

As  indicated  by  Blatt/**  this  "standard  T-dependent  HF  theory"  (which  he  calls  "Husimi  theory") 
is  a  natural  extension  of  zero-temperature  Hartree-Fock  theory  to  finite  temperature. 

However,  it  is  not  the  only  such  extension.  Our  purpose  here  is  to  describe  another,  more 
general,  and  quite  natural  extension.  Recall  that  the  basic  idea  of  zero- temperature  HF  theory 
is  to  look  for  a  single  Slater  determinant  Dq  which  minimizes  the  energy  (D,  HD)/(D,  D)  over  all 
Slater  determinants  D.  Since  the  eigenfunctions  of  the  H  in  Eq.(25)  are  determinants,  it  is 
clear  that  the  minimum  free  energy  computed  with  Eq.  (25)  will,  as  T  -►  0,  approach  Eq  = 

(Dq,  HDq)/(Do,  Dq),  the  correct  zero-temperature  HF  value.  Clearly,  all  that  is  needed  for  this 
is  that  all  the  eigenfunctions  of  H  be  single  Slater  determinants.  But  H,  as  given  in  Eq.  (25),  is 
by  no  means  the  most  general  form  such  that  all  the  eigenfunctions  are  single  determinants.  For 
example,  an  operator  of  the  form  2  e.N.^  +  2  e.  would  also  have  a  complete  set  of  single  de- 

terminantal  eigenfunctions  —  in  fact,  H  would  have  single  determinantal  eigenfunctions  if  it  were 
an  arbitrary  function  E(N^,  .  .  .  )  of  the  occupation  number  operators.  Thus,  we  are  led  to  a 

new  extension  to  finite  temperature  of  zero-temperature  HF  theory  —  choose  the  and  the  func¬ 
tion  E(N^,  N^,  .  .  .  )  to  minimize  the  free  energy  F.  This  theory  is  just  as  natural  as  the  standard 
or  Husimi  extension,  since  it  has  single  Slater  determinants  as  "approximate  energy  eigenfunc¬ 
tions"  ;t  furthermore,  it  can  be  expected  to  be  "better"  in  that  it  will  probably  yield  a  lower  free 
energy.  It  clearly  cannot  give  a  higher  free  energy. 

The  formal  minimization  of  F  in  the  new  context  can  be  carried  out.  It  requires  that  the 
function  E  be 

E(Nf  N2<  .  .  .  )  =  E  (</..,  t*t)  N.  +  }s  v#^) 

-  (^j,  vi/Ki/ij)]  NJNL  (28) 

t  We  feel  that  the  most  general  H  having  its  eigenfunctions  be  Slater  determinants  would  lead  to  a  theory  which  should  be  con¬ 
sidered  the  most  general  extension  to  finite  temperature  of  zero-temperature  HF  theory.  We  have  not  shown  that  the  class  of  IFs 
considered  here  has  that  degree  of  generality. 
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(which  is  just  the  energy  of  the  determinant  defined  by  occupation  numbers  .  .  .  1NL  .  .  .  ),  and  that 
the  ib.  be  solutions  of 


(Ni~Nj)  i ^ 


k*  v^/k>  -  {h*k-  v*k*j>1  (Ni  ~  Nj>  Nk  =  0 


all  i  ^  j 


(29) 


Note  the  similarity  between  these  and  the  standard  HF  equations  [Eq.  (26)];  if  it  were  not  for  the 
terms  in  E  which  are  quadratic  in  the  N's,  the  average  in  Eq.  (29)  (N^  —  N.)  would  factor  into 
(N.  —  N.)  N,  ,  and  cancellation  of  N.  —  N.  would  lead  back  to  Eq.(26)  for  i  =j£j. 

As  we  have  said,  this  theory  completely  corrects  the  physically  serious  deficiency  of  the 
standard  theory  in  the  case  of  infinitely  separated  hydrogen  atoms,  but  we  will  not  go  into  the  de¬ 
tails  here.  At  present,  we  are  studying  this  theory  as  applied  to  hydrogen  atoms  separated  by 


large  but  finite  distances. 


T.  A.  Kaplan 


C.  LASER  SCATTERING  AND  NONLINEAR  EFFECTS 

1.  Hypersonic  Attenuation  in  Quartz  by  Thermal  Brillouin  Scattering 

The  velocity  and  attenuation  of  28-GHz  longitudinal  phonons  along  the  x-axis  of  a -quartz 
have  been  measured  in  a  temperature  range  from  300°  to  600 °K.  High-resolution  Brillouin  scat¬ 
tering  techniques  are  employed  in  this  regime  where  conventional  ultrasonics  techniques  are  not 
applicable.  The  Brillouin  spectrometer  consists  of  a  single-mode  He-Ne  laser  and  a  multiscanned 
confocal  spherical  Fabry-Perot  interferometer.  Photoelectron  count  data  from  a  low  dark  current 
photomultiplier  are  accumulated  in  a  multichannel  analyzer  for  arbitrarily  long  times  to  improve 
the  signal-to-noise  ratio. 

The  velocity,  or  Brillouin  shift,  variation  with  temperature  agrees  within  experimental  er¬ 
ror  with  previous  low-frequency  quartz-controlled  oscillator  measurements.  The  hypersonic  at¬ 
tenuation  is  obtained  from  the  spectral  width  of  the  Brillouin  lines  and  may  be  compared  with 

theoretical  extrapolations  of  ultrasonics  data.  There  is  order  of  magnitude  agreement  with  the 

42  43 

single  relaxation  theories  of  Woodruff  and  Ehrenreich  and  of  Bommel  and  Dransfeld  ;  how¬ 
ever,  both  theories  generally  overestimate  4he  absorption,  and  poorly  represent  the  details  of 
the  temperature  dependence. 

Further  work  is  in  progress  to  extend  the  temperature  range  of  the  measurements. 

A.  S.  Pine 


2.  Light  Scattering  from  Single- Particle  Electron  and  Hole  Excitations 
in  Semiconductors 

Light  scattering  has  been  observed  for  the  first  time  from  screened  single-particle  electron 
excitations  in  GaAs,  InP,  AlSb,  and  CdTe,  as  well  as  from  the  screened  single-particle  hole  ex¬ 
citations  in  p-type  GaAs  doped  with  cadmium,  zinc,  and  magnesium.  The  Raman  scattering  from 

44  45 

plasmons  and  coupled  plasmons  and  phonons,  as  well  as  the  polarization  properties  of  the 
scattered  light  from  these  excitations,  have  been  reported  previously.  Both  these  cases  dealt 
with  the  collective  excitations  of  the  electron  gas  in  a  solid.  Present  results  yield  the  single¬ 
particle  excitation  spectrum  of  the  electron  gas.  The  single- particle  scattering  is  compared  with 

4  6 

effective-mass  calculations  based  on  simple  parabolic  bands,  which  involves  scattering  from 
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electron  density  fluctuations,  as  well  as  a  more  recent  effective-mass  calculation  which  predicts 

47 

a  large  enhancement  due  to  nonparabolicity.  The  experimental  results,  however,  are  in  dis¬ 
agreement  with  both  of  those  theories. 

The  scattered  light  was  generated  using  a  1.06-p  YAG:Nd^+  laser  which  could  operate  with  a 
continuous  output  power  in  excess  of  10  W.  Figure  IV-20  shows  the  Raman  spectrum  of  an  n-type 

GaAs  sample  at  300°,  30°,  and  5°K.  The  peaks  labeled  L  ,  L  and  TO  have  been  described 

45  + 

elsewhere,  and  refer  to  the  coupled  modes  of  the  longitudinal  optical  phonons  and  the  plasmons 

(L  ,  L+),  and  the  transverse  optical  phonon  (TO).  The  scattering  from  the  single-particle  exci¬ 
tation  spectrum  of  the  conduction  electrons  is  seen  at  lower  frequencies  near  the  laser  line.  At 


SHIFT  (cm'h 


Fig.  IV-20.  Raman  spectrum  of  GaAs  (n  =  1.4  X  10^  cm"\  at  300° ,  30°,  and 
5° K.  Scattering  angle  was  90°  with  incident  and  scattered  light  propagating 
along  [100]  crystal  axes. 

T  =  0°K,  the  single-particle  excitation  spectrum  comes  from  electrons  with  momentum  p  that 

are  excited  from  occupied  states  below  the  Fermi  surface  to  unoccupied  states  just  above  the 

Fermi  surface  with  momentum  p  +  q,  where  q  is  the  momentum  change  of  the  scattered  electron. 

46  47 

At  T  =  0°K  and  for  infinite  lifetime,  the  existing  theories  '  predict  a  continuum  spectrum  in¬ 
creasing  linearly  with  w  to  a  peak  near  0.7  qv^,  for  the  parameter  values  involved  here,  and 
cutting  off  at  qvp,  where  v^  is  the  Fermi  velocity.  Samples  with  electron  concentrations  rang¬ 
ing  from  2  X  10«  to  2.8  X  lO1^  cm"3  were  measured.  At  low  temperature,  the  single-particle 
spectrum,  while  smeared  by  a  finite  collision  time,  did  show  a  cutoff  consistent  with  co  =  qv^. 

At  higher  temperatures,  more  electrons  can  participate  in  the  scattering  as  the  electron  distri¬ 
bution  becomes  Maxwellian.  In  the  high-temperature  limit,  the  spectral  shape  of  the  scattered 
47 

light  should  resemble  a  Maxwellian  velocity  distribution,  which  was  found  to  be  consistent  with 
the  experimental  results.  In  the  p-type  GaAs,  scattering  from  both  the  heavy  and  light  hole 
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single-particle  excitations  was  observed  at  300 °K  in  samples  with  hole  concentrations  of  about 
1X10  cm  .  A  more  complete  discussion  of  these  results  is  published  elsewhere.48 

A.  Mooradian 


3.  Steady -State  Thermal  Defocusing  of  Laser  Light 

49 

In  the  previous  Solid  State  Research  Report,  the  time  dependence  of  thermal  defocusing  of 
laser  light  was  discussed  experimentally  and  the  theory  was  presented  80  We  have  carried  out 
further  work  on  the  steady  state  and  the  approach  to  steady  state,  which  includes  a  number  of 
new  experimental  results.  A  steady-state  condition  arises  when  the  temperature  change  due  to 
heat  deposition  is  just  counterbalanced  by  that  due  to  either  heat  conduction  out  of  the  region  of 
the  beam  or  convection  currents  where  material  is  transported  from  the  volume.  In  the  case  of 
convection,  the  hot  material  rises  and  is  replaced  by  cooler  material  from  below.  As  a  result, 
a  vertical  asymmetry  will  develop  for  a  horizontal  laser  beam,  and  the  beam  will  appear  to  be 
bent  downward  and  considerably  distorted.  This  is  to  be  contrasted  with  conduction  which  is 
cylindrically  symmetric,  and  should  therefore  result  in  no  distortions. 

Examining  convection  more  carefully  shows  that  the  acceleration  a,  which  a  particle  expe¬ 
riences  in  the  center  of  the  beam,  will  be  proportional  to  the  total  power  per  unit  volume  ab¬ 
sorbed,  which  is  approximately  a I;  here,  a  is  the  absorption  coefficient,  and  I  is  the  incident 
intensity.  Also,  one  would  expect  convection  to  be  important  when  {i/2.)  X  at^  =  q,  where  t^  is 
the  convection  time  constant,  and  q  is  the  characteristic  vertical  dimension  associated  with  the 
beam.  In  the  case  of  a  horizontal  laser  beam,  q  is  the  beam  diameter  dQ.  From  this  relation, 
it  follows  that 

t^  oc  1  / \TK  oc  i/\f~al  .  (30) 


of  the  time  to  reach  a  maximum  beam  diameter  t  vs  1  /a I  is  presented  for  all  our  CC1.  data. 

max  4 

differs  from  the  convection  time 


The  acceleration  is  also  proportional  to  the  viscosity  of  the  material.  In  Fig.  IV-21,  a  plot 

t  vs  i/a] 
max 

However,  t 

'  max 

by  only  a  constant.  Also  included  is  the  best 
straight  line  with  a  slope  of  l/2,  and  it  is  clear 
that  the  fit  is  quite  good. 

The  theory  for  conduction  effects  was  pre¬ 
sented  by  Gordon,  et  al.,^ *  and  was  discussed 
further  in  Ref.  50  where  it  was  indicated  that  l/t 
should  be  replaced  by  l/t  +  l/t^,  where 


t  = 
c 


d2K 

16  pC 


(31) 


10* 5  10**  10*’ 

1/a  X  (cm3/mW) 

Fig.  IV-21.  A  log-log  plot  of  absorbed  power  per 
unit  volume  al  dependence  of  convection  time,  dem¬ 
onstrating  square-root  dependence  of  t^  on  (al)“"^. 
Liquid  here  is  CCI^,  with  iodine  as  absorber. 


K  =  thermal  conductivity,  p  =  density,  and  = 
heat  capacity  at  constant  pressure.  In  order  to 
study  conduction  effects,  it  is  necessary  to  make 
the  conduction  time  very  much  shorter  than  the 
convection  time.  There  are  at  least  three  ways 
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in  which  this  may  be  accomplished.  The  first  method  is  to  change  the  liquid  from  a  material 
whose  ratio  of  conduction  time  to  convection  time  is  small,  to  one  whose  ratio  is  large.  A  sec¬ 
ond  technique  is  to  change  the  direction  of  propagation  of  the  beam  from  horizontal  to  vertical. 
Now,  instead  of  q  =  dQ,  q  will  be  of  the  order  of  the  cell  length  H,  or  1  /a ,  depending  on  which 
is  shorter.  Thus,  the  convection  time  is  greatly  increased  without  changing  the  conduction  time. 
A  third  method  is  to  use  very  low  absorption  lengths  \/a ,  as  in  Ref.  51.  The  first  two  of  these 
techniques  were  tried,  and  some  results  are  presented  in  Fig.  IV-22.  Data  were  taken  with  a 

vertical  laser  beam  and  a  10-cm  long  vertical  cell  of  both  CS7  (curve  A)  and  CClA  (curve  B).  In 

2  ^  4 

both  cases,  the  intensity  (0.8  W/cm  ),  beam  diameter  (0.294  cm),  and  absorption  length  (9.7  cm) 


SECONDS 

Fig.  IV-22.  A  log-log  plot  of  time  dependence  of  beam  diameter  change  for  10-cm  liquid  cell  containing 
small  amounts  of  iodine,  which  results  in  1/a  =  9.  7  cm.  For  curves  A  (CS2)  and  B  (CCI^,  both  laser  beam 
and  cell  are  vertical.  For  curve  (C),  beam  propagates  horizontally  and  cell  contains  CS^. 

were  the  same.  Also,  in  both  cases,  cylindrically  symmetric  patterns  resulted  for  the  entire 

C  A 

time  development.  Using  the  theory  for  a  parabolic  intensity  profile  gives  the  conduction  times 
of  tc  (CCl^)  =  7.33  sec,  and  tc  (CS^)  =  4.33  sec.  Curve  B  can  be  very  accurately  reproduced 
theoretically  by  using  this  conduction  time  for  CC14  and  replacing  l/t  by  (l/t  +  l/t  ).  Also  given 
in  Fig.  IV-22  is  a  comparison  of  the  vertical  cell  data  for  CS2  (curve  A)  with  that  of  a  horizontal 
cell  with  the  same  constants  (curve  C).  The  exception  is  that  the  input  laser  intensity  was  ~50- 
percent  higher  in  the  horizontal  cell,  causing  the  intercept  to  be  50-percent  larger  in  this  case. 

A  pronounced  ringing  or  underdamped  character  was  displayed  by  the  change  in  beam  diameter 
Ad  for  the  case  of  CS^  in  both  the  vertical  and  horizontal  cell  arrangements.  The  period  of  the 
oscillations  was  much  longer  for  the  vertical  than  for  the  horizontal  cell  and,  in  both  cases,  the 
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oscillations  were  highly  damped.  Apparently,  the  value  of  q  increased  when  switching  from  the 
horizontal  to  vertical  cell,  but  not  by  as  much  as  expected  (l/d ^oi  =  32).  Probably  because  of 
this,  only  the  early  effects  of  conduction  could  adequately  be  predicted  by  the  theory  in  the  case 
of  the  CS^  vertical  cell  data.  More  careful  analysis  of  the  CCl^  data  indicated  that  this  ringing 
was  also  present  in  the  horizontal,  but  not  in  the  vertical,  cell  data.  The  reason  the  theory  does 
not  predict  this  oscillation  undoubtedly  lies  with  two  approximations.  One  is  that 


f 


I(tf)  dt'  =  It 


(32) 


If  the  past  history  of  the  intensity  profile  is  included  with  the  instantaneous  variations,  an  effec¬ 
tive  time  lag  should  result,  thus  causing  the  oscillations.  A  second  reason  is  that  replacing  l/t 
by  (l/t  +  l/t  )  is  correct  only  to  first  order,  and  higher  order  terms  may  be  significant  for  some 

situations.  „  T  ^ 

R.  L.  Carman 

P.  L.  Kelley 


4.  Single-Mode  Operation  of  a  High-Power  Pulsed  Ruby  Laser 

a.  Transverse-Mode  Selection 

Our  most  reliable  transverse-mode  selection  technique  employed  a  diamond  pinhole  and  two 
long  focal  length  lenses,  as  illustrated  in  Fig.  IV-23(a).  The  resulting  cavity  length  is  approxi¬ 
mately  1.5  meters.  The  output  beam  from  this  configuration,  using  the  prism  as  a  Q-switch,  is 


Fig.  IV-23.  (a)  Schematic  illustrating  experimental  arrangement  used  for  transverse¬ 
mode  selection,  (b)  Interferogram  of  non-Q-switched  longitudinal -mode  structure 
made  with  a  Fabry-Perot  etalon  (inter-order  spacing  11.8GHz,  resolution  ~600  MHz). 


a  4-mm-diameter  beam  with  a  0.5-  (theoretical  limit)  to  1-mrad  divergence,  and  a  power  of  0.1 
to  1  MW  which  is  limited  by  air  breakdown  in  the  pinhole.  When  the  prism  was  fixed  and  Q- 
switching  was  accomplished  with  a  passive  dye,  then  peculiar  near-field  patterns  as  well  as  air 
breakdown  in  the  pinhole  were  unavoidable,  even  with  total  output  power  as  low  as  50 kW.  The 
air  breakdown  appears  to  produce  a  plasma  which  is  opaque,  and  thus  cuts  off  further  lasing. 

In  the  case  of  the  passive  Q-switch,  a  second  method  of  transverse-mode  selection  was  also 
tried;  this  consisted  of  removing  the  pinhole  and  placing  a  second  aperture  (optional  aperture)  as 
indicated  in  Fig.  IV-23(a).  Both  apertures  were  identical  in  size  and  centered  on  a  common  axis. 
Using  4-mm-diameter  apertures,  the  laser  output  consisted  of  three  beams,  each  approximately 
1  mm  in  diameter  with  about  1 . 5-mrad  divergence.  It  was  concluded  that  the  passive  dye  used  in 
conjunction  with  a  flame  fusion-grown  ruby  produces  spatial  inhomogeneities  whose  characteristic 
dimensions  are  too  small  to  render  these  beams  useful. 
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For  the  case  of  the  pinhole,  the  existence  of  several  diffraction-limited  hot  spots  could 
result  in  power  densities  high  enough  to  cause  localized  air  breakdown  and,  from  the  ensuing 
scattering  and  absorption,  a  peculiar  output  near-field  pattern.  Using  the  passive  dye  as  a 
mode-locking  element,  the  requirements  on  spatial  inhomogeneities  may  well  be  too  high  for  re¬ 
liable  operation  of  a  diffraction-limited  mode-locked  laser. 


b.  Longitudinal-Mode  Selection 


For  non-Q-switched  operation,  the  arrangement  in  Fig.  IV-23(a)  produced  sufficient  mode 
selection  near  threshold  to  allow  for  operation  with  linewidths  of  less  than  600  MHz,  the  resolu¬ 
tion  limit  of  our  present  Fabry-Perot  interferometer;  i.e.,  less  than  six  modes  were  present. 

A  typical  interferogram  is  shown  in  Fig.  IV-23(b).  An  interferometer  with  a  resolution  less  than 
100  MHz  is  now  being  constructed  to  confirm  single-mode  operation. 

When  the  laser  was  Q-switched  by  rotating  a  prism,  the  scheme  illustrated  in  Fig.  IV-23(a) 
was  insufficient  for  producing  single-longitudinal-mode  operation.  By  adding  a  second  feedback 
arrangement,  which  comprised  the  second  arm  of  an  inverted  Michelson  interferometer,  single¬ 
mode  operation  was  obtained.  The  45°  beamsplitter  had  a  60-percent  reflectivity,  while  the  re¬ 
maining  mirror  had  a  100-percent  reflectivity.  The  only  mode  which  will  oscillate  in  this  scheme 
is  the  one  for  which  reflection  from  the  100-percent  mirror  destructively  interferes  with  the  light 
directly  reflected  out  of  the  cavity  by  the  beamsplitter.  Reliable  longitudinal-mode  selection  re¬ 
quires  stable  mirror  mounts;  we  are  now  investigating  this  problem. 

R.  L.  Carman 
F.  H.  Perry 
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